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Abstract 

The transformation of a Laguerre series f(z) = YJ^= ^ (^) to ^ power series 

f( z ) = L^=0 7' !Z " i s discussed. Since many nonanalytic functions can be expanded in 
terms of generalized Laguerre polynomials, success is not guaranteed and such a transfor- 
mation can easily lead to a mathematically meaningless expansion containing power series 
coefficients that are infinite in magnitude. Simple sufficient conditions based on the decay 
rates and sign patters of the Laguerre series coefficients A.j°^ as n — > °o can be formu- 
lated which guarantee that the resulting power series represents an analytic function. The 
transformation produces a mathematically meaningful result if the coefficients either 

decay exponentially or factorially as n — » °o. The situation is much more complicated - 

(ft) 

but also much more interesting - if the A„ decay only algebraically as n — > °°. If the 

(ft) 

ultimately have the same sign, the series expansions for the power series coefficients 

(ft) 

diverge, and the corresponding function is not analytic at the origin. If the X„ ultimately 
have strictly alternating signs, the series expansions for the power series coefficients still 
diverge, but are summable to something finite, and the resulting power series represents an 

analytic function. If algebraically decaying and ultimately alternating Laguerre series co- 

(ct) 

efficients A,, possess sufficiently simple explicit analytical expressions, the summation 
of the divergent series for the power series coefficients can often be accomplished with the 

(ft) 

help of analytic continuation formulas for hypergeometric series p+\F p , but if the A„ 
have a complicated structure or if only their numerical values are available, numerical 
summation techniques have to be employed. It is shown that certain nonlinear sequence 
transformations - in particular the so-called delta transformation [E. J. Weniger, Comput. 
Phys. Rep. 10, 189 - 371 (1989), Eq. (8.4-4)] - are able to sum the divergent series oc- 
curring in this context effectively. As a physical application of the results of this article, 
the legitimacy of the rearrangement of certain one-range addition theorems for Slater-type 
functions [I. I. Guseinov, Phy. Rev. A 22, 369 - 371 (1980); Int. J. Quantum Chem. 81, 
126- 129 (2001); Int. J. Quantum Chem. 90, 1 14 - 1 18 (2002)] is investigated. 

PACS numbers: 02.30.Gp, 02.30.Lt, 02.30.Mv, 02.60.-x 

AMS classification scheme numbers: 30B10, 30B40, 33C45, 40A05, 65B10 
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1 Introduction 

ice) 

The generalized Laguerre polynomials L„ (z), whose most relevant properties are reviewed 
in Section |2] are a very important class of orthogonal polynomials with numerous mathemat- 
ical and physical applications. There is an extensive literature both on their mathematical 
properties as well as on their applications, and any attempt of providing a reasonably com- 
plete bibliography would be hopeless. Let me just mention that the radial parts of bound state 
hydrogen eigenfunctions and of several other physically relevant complete and orthonormal 
sets of function / : M 3 — » C are essentially generalized Laguerre polynomials (see for example 
11 1 351 Sections IV and V] and references therein). 

As discussed in more detail in Section [2] the generalized Laguerre polynomials (z) 
form a complete orthogonal polynomial system in the weighted Hilbert space L 2 ag _ z ([0,°°)) 
defined by (12.7b . which is based on an integration over the interval [0,°°) involving the weight 
function w(z) = z"exp(— z). Accordingly, the topic of this article are infinite Laguerre series 
of the following type: 

f(z) = £ A„ (a) L< a) (z), (1.1 a) 

n=0 

X " a) = rtnt' In r^e-4 a) (z)/(z)dz. (Lib) 
i(a + n + 1) Jo 

The expansion coefficients A„ are essentially inner products, utilizing the orthogonality and 
the completeness of the generalized Laguerre polynomials in £ 2 a e -z ( [0, <*>)). 

Not all expansions in terms of generalized Laguerre polynomials can be interpreted in 
a Hilbert space setting. There are expansions of special functions in terms of generalized 
Laguerre polynomials with variable and index-dependent superscripts. A simple example is 
the generating function 11021 p. 242] 

£ L [ ,r n \x)t n = e-*(l +0", |*| < 1. (1-2) 

n=0 

Expansions of that kind cannot be derived via a straightforward application of the orthogonal- 
ity of the generalized Laguerre polynomials. In the article by Sanchez-Ruiz, Lopez- Artez, and 
Dehesa II 131 . transformation formulas for generalized hypergeometric series were employed 
for the construction of Laguerre expansions with variable superscripts. However, nonorthog- 
onal Laguerre expansions of that kind are not the topic of this article. 

It is generally accepted that orthogonal expansions are extremely useful mathematical 
tools and that they have many highly advantageous features. This is, however, not the whole 
truth, in particular if we want to approximate functions. There are situations, in which alter- 
native representations are (much) more useful. Hilbert space theory only guarantees that an 
orthogonal expansion converges in the mean with respect to the corresponding norm, but not 
necessarily pointwise or even uniformly. Therefore, orthogonal expansions are not necessarily 
the best choice if we are interested in local properties of functions. 

Convergence in the mean is - although completely satisfactory for many purposes such 
as the evaluation of integrals - a comparatively weak form of convergence. In practice, it is 
therefore often desirable or even necessary to construct alternative representations possessing 
more convenient properties (see also the discussion in [46 1). 

A very desirable feature of functions /: C — * C is analyticity in the sense of complex 
analysis. This means that a function / can be represented in a neighborhood of the origin by 
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a convergent power series, 

/(z) = X>z", (1.3) 

n=0 

and that the coefficients y„ of this power series essentially correspond to the derivatives of / 
at z = 0. 

If the function / defined by the Laguerre series (II . lb is explicitly known, it is usually not 
too difficult to decide whether / is analytic or not, and if / is analytic it is normally not too 
difficult to construct at least the leading terms of its power series expansion (11.31 l. 

Unfortunately, the situation is not always so good, and it can even happen that only the 
numerical values of a finite number of expansion coefficients xjf in (II. Il l are known. In 
such a case, it would certainly be helpful if we could relate the properties of the expansion 
coefficients X„ - in particular their decay rate and their sign pattern - to the analyticity of 
the function f(z) defined by the Laguerre series (II. lb . It would also be helpful if we could 
construct at least the leading power series coefficients y„ - either exactly or in an approximate 
sense - from the coefficients A„ of the Laguerre series ( II. 11 1. . 

Pollard 11071 . Szasz and Yeardley [125|, and Rusev 111 1 11 had investigated Laguerre ex- 
pansions of the type of ( II. Il l of analytic functions and analyzed their regions of convergence. 
However, the inverse problem - the formulation of criteria which guarantee that a Laguerre 
expansion of the type of (II. Il l represents a function analytic in a neighborhood of the origin - 
seems to be essentially unexplored. I am only aware of short remarks by Gottlieb and Orszag 
||49l p. 42] and by Doha ATI p. 5452], respectively, who stated without detailed proof that 
such a Laguerre series converges faster than algebraically if the function under consideration 
is analytic at the origin. 

Many Laguerre series are known which seem to confirm the claim of Gottlieb and Orszag 
and of Doha, respectively. The probably most simple example is the well known generating 
function (1 — t )~ a_1 exp(zf /[t — 1]) of the Laguerre polynomials. For f | < 1, it is an entire 
function, and the coefficients of its Laguerre series (I5.81 l decay exponentially for \t \ < 1. 

Another example is the Laguerre series (I3.31 l for the general power function z p with nonin- 
tegral p G M.\No- Obviously, z p with p e R\No is not analytic at z = 0. The series coefficients 
in ( 13.31 1 decay algebraically and possess for sufficiently large indices the same sign. As shown 
in Section [3] it is nevertheless possible to transform the Laguerre series d3.31 > for z p with 
p G K \ No to the formal power series ( 13.7b . However, the formal power series (13.71) is not a 
mathematically meaningful object: For sufficiently large indices, the power series coefficients 
in (13.71 ) are according to ( 13.9b all infinite in magnitude. Thus, the Laguerre series (13.3b for z p 
with p €R\Nq also agrees with the claim of Gottlieb and Orszag and of Doha, respectively. 

If, however, we augment the Laguerre series ( 13.3b for z p by an additional alternating sign 
(—1 )", we obtain the Laguerre series ( 14.7b . Its transformation to a power series yields accord- 
ing to ( 14.12b a confluent hypergeometric function i Fi ( — p ; a + 1 ; z/2) , which for Re ( a ) > — 1 
is an entire function. The coefficients in the Laguerre series ( 14.7b decay algebraically in mag- 
nitude. This implies that the transformation of ( 14.7b to a power series leads to power series 
coefficients y„ that are represented by divergent series expansions. However, the alternating 
sign (—1)" in ( 14.7b makes it possible to associate finite values to these divergent series, i.e., 
the series expansions for the power series coefficients are now summable. Thus, we obtain a 
function, which is analytic at z — and which does not agree with the claim of Gottlieb and 
Orszag and of Doha, respectively. 

The example of the closely related Laguerre series ( 13.3b and ( 14.7b shows that it makes a 

(a) 

huge difference if Laguerre series coefficients A« , that decay algebraically in magnitude, 
ultimately have strictly alternating or strictly monotone signs. Thus, the claim of Gottlieb and 
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Orszag and of Doha, respectively, is imprecise and ignores the pivotal role of divergent, but 
summable series expansions for the power series coefficients y„ in (11.31 ). Some basic facts 
about the summation of divergent series are reviewed in AppendixlAl 

A detailed investigation of the analyticity of functions represented by Laguerre expansions 
of the type of (11.1b is the topic of this article. The central result is the transformation formula 
(13. 14b , which yields a formal power series expansion for a function represented by a Laguerre 
expansion. The crucial question is whether the inner fi series in ( 13.141 ) for the power series 
coefficients y„ converge. If these fi series do not converge and are also not summable to 
something finite, then the function under consideration is not analytic at the origin and the 
resulting formal power series that is mathematically meaningless since contains coefficients y„ 
which are infinite in magnitude. Consequently, it is comparatively easy to relate the analyticity 

of a function / represented by a Laguerre series at the origin z = with the decay rate and the 

( ex) 

sign pattern of its coefficients A„ . 

Analytical manipulations can only lead to closed form expressions for the power series 
coefficients y„ if the Laguerre series coefficients jM possess a sufficiently simple structure. 

This undeniable fact may create the false impression that the formalism developed in this 

( tx) 

article is restricted to Laguerre series with very simple coefficients A„ , and that this for- 
malism is at best suited for an alternative rederivation of known generating functions of the 
generalized Laguerre polynomials. 

As discussed in Section [6j it is, however, often possible to construct numerical approxi- 
mations to the leading power series coefficients y„ from the coefficients X„ of the Laguerre 
series, even if the inner fi series in (13.14b diverges. The necessary summations can be done 
effectively with the help of certain nonlinear sequence transformations - Wynn's epsilon al- 
gorithm ( IB. 4b and the two Levin-type transformations (IB. 14b and ( IB. 15b - whose basic prop- 
erties are reviewed in Appendix 151 While Wynn's epsilon algorithm ( IB. 4b . which produces 
Pade approximants if the input data are the partial sums of a power series, is now fairly well 
known among (applied) mathematicians and theoretical physicists, the Levin-type transfor- 
mations ( IB. 14b and (IB. 15b . which in some cases were found to be remarkably powerful, are 
not nearly as well known as they deserve to be. 

As a physical application of the mathematical formalism developed in this article, the 
legitimacy of rearrangements of certain one-range addition theorems for Slater-type functions 
with in general nonintegral principal quantum numbers performed by Guseinov ll53l l54l l56l 
is investigated in Section [7] It is shown that the one-center limits of Guseinov's rearranged 
addition theorems do not exist if the principal quantum numbers of the Slater-type functions 
are nonintegral. 



2 Generalized Laguerre Polynomials 

The generalized Laguerre polynomials L„ (z) with n £ No and Re(oc) > — 1 are orthogonal 
polynomials associated with the integration interval [0,°°) and the weight function w(z) — 
z"exp(-z). 

In this article, the mathematical notation for generalized Laguerre polynomials (see for 
example [ 102, Chapter 5.5]) is used. Additional conventions, as used predominantly in older 
books and articles on quantum theory, were discussed by Kaijser and Smith [87, Footnote 1 
on p. 48] and also in 11 1 5 II Section 4]. 

The generalized Laguerre polynomials possess an explicit expression as a terminating 
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confluent hypergeometric series HI 021 p. 240] 



Ll a \z) = { -^±^ lFl (-n-a + l;z), (2.1) 



and they can also be defined via their Rodrigues relationship Ml 021 p. 241] 

In the special case a = 0, it is common to drop the superscript since we obtain the usual 
Laguerre polynomials 111021 p. 239]: 

L„(z) = 4 0) (z), neN Q . (2.3) 

Either the explicit expression (12. U or the Rodrigues relationship ( 12.2b can be used to define 
generalized Laguerre polynomials with essentially arbitrary complex values of the superscript 
a. But in the orthogonality relationship of the generalized Laguerre polynomials II 1021 p. 
241], 

T(a + n + l) 



z a e- z LX>(z)Lf>(z)dz = y — ' 8, nn , m,»eN , (2.4) 

n! 

it is necessary to assume Re(a) > — 1 because otherwise this integral does not exist. However, 
as discussed by Kochneff [95 1, this condition can be relaxed if this and related integrals are 
reinterpreted as Hadamard finite part integrals. 

In the vast majority of calculations involving generalized Laguerre polynomials, the su- 
perscript a is real and also positive. In practice, more general complex values of a play a 
negligible rule. Therefore, in the following text it is always tacitly assumed that all integrals 
exist in the ordinary sense, which avoids the complications of Hadamard finite part integrals, 
and that the the superscript a of a generalized Laguerre polynomial is a real number satisfying 
a > -1. 

The orthogonality relationship ( 12.4b suggests the introduction of the following inner prod- 
uct containing the weight function w(z) = z"exp(— z) for functions f,g: C — > C: 



(f\g) z « e -*,2 = L z a ^[f(zWg(z)dz. (2.5) 
This inner product gives rise to the norm 



ll/IU-2 - V(/l/) Z «e-*,2- (2 ' 6) 

We then obtain the following Hilbert space of square integrable functions: 

4<e-;([0H) - {/: C^C| ^°Ve--'|/(z)| 2 dz<-} 

= {/: C^C|||/|| z « e -z, 2 <oo}. (2.7) 

The completeness of the generalized Laguerre polynomials in this weighted Hilbert space 
is a classic result of mathematical analysis (see for example ll79l p. 33], Ml 141 pp. 349 - 351], 
or [1 1 301 pp. 235 - 238]). Thus, the normalized polynomials 



T(a+n + \) 



1/2 

zi°°(z), »gN , (2.8) 
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are a complete and orthonormal polynomial system in L?„ e _ z ([0,°°)), satisfying 

z a e~ z J?L a) (z) (z) dz = 8 m „ . (2.9) 

Accordingly, every / 6 l? a _ z ([0,°°)) can be expanded in terms of the normalized polynomi- 

als{^(z)}^ : 

f{z) = £% (o) JSf» (o) (z), (2.10a) 

n=0 



,(a) 



z a e- z ifr ) (z)/(z)dz. (2.10b) 



This expansion in terms of the normaized polynomials Jtf„ ' (z), which converges in the mean 
with respect to the norm (12.6b . is nothing but the Laguerre series ( II. lb in disguise. In agree- 
ment with ( 12.8b . we only have to choose 



r(a + « + l) 



1/2 



A„ (a) (2.11) 



to see that the expansions (II . lb and ( 12.10b are identical. 

Let me emphasize once more that Laguerre expansions converge in general only in the 
mean, but not necessarily pointwise (see for example [2]). Additional conditions, which a 
function has to satisfy in order to guarantee that its Laguerre expansion also converges point- 
wise, were discussed by Szego II 1261 Theorem 9.1.5 on p. 246]. 

Hilbert space theory can be used to derive sufficient criteria, which the coefficients A" 
or have to satisfy in order to guarantee that the equivalent expansions ( II. lb and (12.10b 
converge in the mean. The basic requirement is that the norm ( 12.6b of both / and its Laguerre 
expansion (12.10b must be finite: 

ll/IU-2 = (L|^ (a) | 2 ) 1/2 < °°- (2-12) 
It follows from ( 12.1 lb that this condition can also be reformulated as follows: 

ina±^ m y\„ (213) 

«=o n - > 



\z a e- z ,2 



A sufficient condition, which guarantees that the infinite series LJJLol"^" ^ | 2 m (12.12b 
converges, is 

\%l a) \ 2 ~ n- l - £ , n^oo, e >0. (2.14) 

It follows from d2.1 lb and (12.14b that the Laguerre series ( 11.11 ) converges in the mean if its 
coefficients Xj,^ satisfy 

r(a + rc+l)| («)|2 _i_ £ „ 

— - An ; ~ n n^°°, e>0. (2.15) 

n! 1 1 

Now, we need the asymptotic approximation [1 , Eq. (6.1.47) on p. 257] 

r(z + a)/r{z + b) = z a - h [1+0(1/ z)], z^~, (2.16) 
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which is the leading term of an asymptotic expansion first derived by Tricomi and Erdelyi 
11 1 3 1 II that holds - as emphasized by Olver [104, p. 119] - without restrictions on a,b G C. 
Thus, we obtain: 

r(q + « + l) ^ n q +o(w q- 1)i (21?) 

i (cc) I 

This translates to the sufficient convergence condition that | X„ | must decay like 

|Ai a) | ~ «-[ a+£+1 ]/ 2 , «^oo, e >0, (2.18) 

or faster. If this asymptotic condition is satisfied, then the Laguerre series ( II. Il l for f(z) 
converges in the mean with respect to the norm ( 12.6b of the Hilbert space L^- : ^ a ([0,°°)) . 

The sufficient convergence condition ( 12.181 ) shows that the coefficients A„ can decay 
extremely slowly, which translates to a possibly extremely slow convergence of Laguerre 
series of the type of (II. It . Thus, it would be overly optimistic to assume that all Laguerre 
series are necessarily computationally useful. 



3 The Laguerre Series for the General Power Function 

In the theory of orthogonal expansions, which converge in the mean with respect to the norm 
|| • || of some Hilbert space the decisive criterion is that the function, which is to be 
expanded, has to belong to . Consequently, any function /: C — > C with finite norm 
(I2.61 l belongs to the weighted Hilbert space L 2 a _ z ([0,°°)) and can be expanded in terms of 
generalized Laguerre polynomials according to ( 11.11 ). The resulting expansion ( II. It converges 
in the mean with respect to the norm (12.61 ) of the Hilbert space Z/?a e - : ([0,°°))- 

The existence and convergence of a Laguerre series of the type of dl.lt does not allow 
any conclusion about the analyticity of the corresponding function. The weighted Hilbert 
space L? a ; ([0,oo)) contains many functions that are obviously not analytic at the origin. 
Therefore, attempts of constructing a power series from the Laguerre series can easily lead to 
a disaster. The complications, which can occur in this context, can be demonstrated via the 
general power function z p with nonintegral p G R \ No- 

In integrals over the positive real semi-axis, the weight function z a e~ z becomes at least 
for Re(a) > very small both as z — > and as z — > °°. Consequently, z a e~ z suppresses 
in the inner product (12.5b the contribution of the remaining integrand for small and large 
arguments. Therefore, the general power function z p with p G R\Mq possesses at least for 
sufficiently large values of a an expansion of the type of dl.lt in terms of generalized Laguerre 
polynomials. The existence of this expansion is guaranteed if z p belongs to the Hilbert space 
L?a e - Z (|0,°o)), or equivalently, if 



j; 



7 a+2p 



< °o (3.1) 



holds. Thus, we have to require that a + 2p > —1 holds. 

For the construction of a Laguerre series for z p , we only need the integral 11501 Eq. 
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(7.414.7) on p. 850] 

'e-V4 a) (f)df 

■s-^- l z F\{-n^ + V,a+\;\/s) 



r(/3 + l)r(a + « + l) __ p _ K 



n!r(a + l) 

Re(jB)>-l, Re(i)>0, (3.2) 



to obtain after some essentially straightforward algebra [43, Eq. (16) on p. 214]: 

S^-W^i^^ »**\">. «+*»-•■ <«> 

If we set p = m with m G No, then the infinite series on the right-hand side terminates because 
of the Pochhammer symbol (— m) n and we obtain the following finite sum (see for example 
lfT09lEa. (2) on p. 207]): 

z m = (a + 1),„ £ L (°0( z) ( m e N(J) a + 2 m > -1 . (3.4) 

„=o (« + !;« 

Although intimately related, there are nevertheless some fundamental differences between 
the two Laguerre expansions d33T > and (EUt . The finite sum formula ( I3.4l i is a relationship 
among polynomials and therefore certainly valid pointwise for arbitrary z G C as well as 
analytic at the origin in the sense of complex analysis. 

In the case of the infinite series expansion ( 13.3b . we only know that it converges in the 
mean with respect to the norm d2.61 >. but we have no a priori reason to assume that this 
expansion might converge pointwise for arbitrary p G R\No- Moreover, z p with p G R \ No 
is not analytic at the origin. 

The validity of the finite sum ( 13.4b can be checked by inserting the explicit expression 
( 12.11 ) for the generalized Laguerre polynomials. If we then rearrange the order of the two 
nested sums, we obtain after some algebra: 

Next, we use the relationship [ 1 , Eq. (3.1.7) on p. 10]: 

E = «*. »eN , (3.6) 

for binomial sums, which shows that the inner sum in ( 13.5b vanishes unless we have k = m. 
Thus, we only need (— l) m (— m) m = (l) m = ml to arrive at the trivial identity z m = z m which 
proves the correctness of (13.4b . 

In the case of the infinite series ( 13.3b for z p with p G R \ No, we can also insert the explicit 
expression ( 12. lb for the generalized Laguerre polynomials into it and rearrange the order of 
summations. We then obtain after some algebra: 

ri^ + i) t=o [a + l) k k\ 



Ernst Joachim Weniger: On the Analyticity of Laguerre Series 



3: The Laguerre Series for the General Power Function 



10 



Superficially, it looks as if we succeeded in constructing a power series expansion for the in 
general nonintegral power z p . However, the generalized hypergeometric series \Fq with unit 
argument is the limiting case z — ► 1 of the so-called binomial series HI 021 p. 38]: 

lFo (a;z) = t (~ a )(-z) m = (l-zT a , kl<l- (3-8) 

m=0 \ m / 

If we set a = k — p with k € No and p £ M \ No, we obtain for the hypergeometric series iFo 
inlEzJ: 

P < 0, 

0, k<p>0, (3.9) 
oo, k > p > 0. 

Thus, disaster struck and the power series ( 13.7b contains infinitely many series coefficients 
that are infinite in magnitude. 

I am aware of an article by Villani 013311 who tried to make sense of perturbation expan- 
sions with divergent terms. Since, however, this is the only article on this topic, which I am 
aware of, I am tempted to believe that Villani's approach was not overly fertile. Therefore, I 
will stick to the usual mathematical convention that a power series with coefficients, that are 
infinite in magnitude, does not exist as a mathematically meaningful object. 

It is important to note that series with divergent terms and divergent series are not the 
same. In the case of divergent series, all terms are finite, but the conventional process of 
adding up the terms successively does not lead to a convergent result. Nevertheless, it is often 
possible to associate a finite value to a divergent series with the help of a suitable summation 
technique. As reviewed in Appendix|A] divergent series and their summation have been and 
to some extent still are a highly controversial topic. The summation of divergent series plays 
a major role in Sections (0]l and ©. 

Since the hypergeometric series \Fo(k — P\ 1) m (13.7b does not converge for all indices k, 
the Laguerre series (13.3b for z p with p € R \ No cannot be reformulated as a power series in 
z by interchanging the order of summations. Of course, this is a mathematical necessity: The 
general power function z p with p £ R \ No is not analytic at the origin, which implies that a 
power series about z = cannot exist. 

So far, the analysis of this Section has only produced obvious results and no new insight: 
The integral power z m with m £ No is analytic at the origin. Consequently, it must be possible 
to reformulate its finite Laguerre expansion (13.4b as a polynomial in z. In contrast, the non- 
integral power z p with p £ M. \ No is not analytic at the origin. Accordingly, a power series 
for z p about z = cannot exist. At least formally, the infinite Laguerre expansion (13.3b for 
z p can be rearranged to yield the power series ( 13.7b . but this power series is mathematically 
meaningless since almost all of its series coefficients are infinite in magnitude. 

Nevertheless, the example of the general power function z p is instructive since it shows 
that mathematics cannot be cheated by formally rearranging Laguerre series. Therefore, we 
can try to use the approach described above also in the case of essentially arbitrary infinite 
Laguerre series of the type of ( II. lb . We insert the explicit expression ( 12.1b of the generalized 
Laguerre polynomial into the Laguerre series for some function, and rearrange the order of 
summations. As a final step, we have to analyze whether and under which conditions the inner 
infinite series expansions for the coefficients of the resulting the power series converge. 

But first, let us consider a partial sum of the general Laguerre series dl.lt : 

f N (z) = E ti a) Ll a \z), N£N . (3.10) 

n=0 
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Such a finite sum is simply a polynomial in z and it is always possible to reformulate it by in- 
terchanging the order of the nested finite summations. If we insert the explicit expression ( 12.1b 
of the generalized Laguerre polynomials into (13.10b and rearrange the order of summations, 
we obtain after some algebra: 

f M _ f a W (« + !)» y (~")v ^ mn 

Mz) -„ E o A " ^^ v E (^TTj;^ (311) 
= E^nE ^±I^>. (3 , 2) 

v=o (a + l)vV! „f v «! 
This expression can be streamlined further, yielding 

fM-Y { ~ z)V V ( g + v + 1 )M , l W n IT, 

MzJ - L —^7- L a m+v ( 3 - 13 ) 

v=0 v - /i=0 ^' 

Thus, in the case of finite Laguerre expansions of the type of ( 13.10b . a rearrangement of the 
order of the nested finite summations is always possible. 

Let us now consider the rearrangement of the infinite series ( 11.1b . By inserting the explicit 
expression ( 12. lb of the generalized Laguerre polynomials into dl.lb and rearranging the order 
of summations, we formally obtain the following power series in z: 

m _ £ (^T £ l«±^k i&t . (3 . 14) 

v=0 v - fi=0 

If we compare (I3.13I , and (13.14b . we immediately see that the rearrangement of an infinite 
Laguerre expansion is not necessarily possible since we now have an inner infinite series 
instead of an inner finite sum. Accordingly, many things can go wrong if we mechanically 
perform the limit N — > °° in ( 13.13b . The power series (13.14b for f(z) makes sense if and only 
if the inner series on the right-hand side of ( 13.14b converges for every v G No- Otherwise, we 
have a formal power series with expansion coefficients that are infinite in magnitude. This 
scenario corresponds to the formal, but mathematically meaningless power series (13.7b for z p 
wifhp GK\N . 



4 Algebraically Decaying Series Coefficients 

In this Section, the convergence of the inner /I series in the rearranged Laguerre expansion 
( 13.14b is analyzed by making several assumptions on the large index (n — » 00) asymptotics of 
the coefficients k^ a \ 

The sufficient convergence condition ( 12.181 ) shows that the coefficients A„ in ( II. lb can 
decay algebraically in n, which in practice implies (very) bad convergence. Thus, let us 
assume for the moment that the X„ all have the same sign at least for sufficiently large 
indices n, and that they possess the following large index asymptotics: 

X,{ a) ~ n-P, n^oo, p >o. (4.1) 

For an analysis of the convergence of the inner fi series, it is helpful to rewrite ( 13.14b as 
follows: 

fM _ 1 y (--*) V f. r(q + /x + v + l) (a) 

f(x) r(a+I) h Ja + lj^. £i M+v ■ (4 ' 2) 
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We first analyze the large index asymptotics of the factor F(a + ji + v + 1)/ With the help 
of (12.161 1. we obtain for fixed and finite v € No the following asymptotic approximation: 

r(«+ M +v+i) ^ a+V] ^ (4J) 

For fixed and finite v G No, the asymptotic estimate ( 14.11 ) translates to 

4+v ~ (M + V)-" = /i-" + O^-^ 1 ) , (4.4) 
Combination of (14.3b and ( I4.4l i yields 
r(a + jii + v + l) 



A^ v ~M a+v - p , M-°°- (4-5) 



Ml 



Thus, the inner /I series in ( 13.14b diverges at least for sufficiently large values of the outer 
index v if the series coefficients A« occurring in (13.14b ultimately have the same sign and 
decay algebraically like a fixed power j3 of the index n. Thus, a function represented by a 
Laguerre expansion with ultimately monotone and algebraically decaying series coefficients 
cannot be analytic in a neighborhood of the origin. This conclusion is in agreement with the 
remarks by Gottlieb and Orszag ll49l p. 42] and by Doha iPTO p. 5452], respectively, who had 
stated that the Laguerre series for a given function converges faster than algebraically if the 
function under consideration is analytic at the origin. 

As discussed in more details in Section [5] the general power function z p with p G K \ No 
is not analytic at the origin. This fact can also be deduced from its Laguerre series ( 13.3b . 
With the help of (12.16b . we obtain the following leading order asymptotic estimate for the 
coefficients in (13.3b : 

r(-p + n) -a-p-t n ^ QO (46) 

T(a + « + l) 

Comparison with (12.18b shows that this asymptotic estimate implies the convergence of the 
Laguerre series (13.3b forz p with respect to the norm (12.6b of the Hilbert space I? ae ^ x ([0,°°)) 
if a + 2p > — 1 holds. However, this estimate also shows that z p with p £ K \ No cannot be 
analytic at the origin z = 0. 

The convergence properties of monotone and alternating series differ substantially. A 
monotone series Y^=o a n, whose terms all have the same sign, converges, if the series terms 
a n decay at least like a n — 0{n~'~ £ ) with e > as n — > °° or faster. In contrast, an alternat- 
ing series £" = o(— \)"\b n \ converges if the terms b„ decrease in magnitude and approach zero 
as n — > oo. This alone would not suffice to guarantee the convergence of the inner fi series 
in ( 13.141 ). However, alternating series have the undeniable advantage that summability tech- 
niques for divergent series can be employed (further details as well as numerous references 
can be found in Appendices|A]and|B]). In this way, it is frequently possible to associate a finite 
value to a divergent alternating series 5Z.JT=o ( — 1)"|^b|> whose terms do not vanish as n — > °° 
and even grow in magnitude with increasing index. 

Thus, it should be interesting to investigate whether Laguerre series of the type of ( II. lb 

f eel 

with strictly alternating series coefficients A„ , whose absolute values are at least asymptot- 
ically proportional to a fixed power of the index n, correspond to functions that are analytic in 
a neighborhood of the origin. 

As an example of a Laguerre series with ultimately strictly alternating and algebraically 
decaying coefficients, let us consider the following Laguerre series, which differs from the 
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Laguerre series (13.31 ) for the general power function z p only by the alternating sign (— 1 )": 

(a) _ T(p + a + l) ~ (-p)„ (g) 

Gp (Z) - r(a + l ,^ ( 1} (a + l)> 

peM\N , a+2p>-l. (4.7) 

The convergence condition a + 2p > — 1 guarantees that G p (z) belongs just like the power 
function z p with p 6 IR\No to the Hilbert space L 2 a _. v ([0,°°)). The case p G No is excluded 

fee) 

because then G p would be a polynomial in z, whose analyticity at the origin is obvious. 

(ot) 

It is, however, unclear whether G p (z) with p G M \ No is analytic at the origin, i.e., 
whether the inner fi series in the rearranged Laguerre series (13.141) converges for arbitrary 
values of v £ No. For an investigation of this question, let us define 

tf) = ( _ iy .r(p + «+i) n-p+n) 
y ' r(-p) r(«+« + i) 

Inserting this into the modified rearranged Laguerre expansion ( 14.21 ) yields: 

The generalized hypergeometric series \Fq in (14.10b is a special case of the binomial series 
( T3~8l) which can be expressed in close form. Thus, we obtain: 

iF (-p + V;-l) = Urn (l-z) p ~ v = 2 P ~ V . (4.11) 

Since ( 12.16b implies (— p + v)^/ju! = 0(/X v ~''~ 1 ) as ji — > °o, the hypergeometric series 
ii*b(— p + v;z) converges only for |z| < 1. Thus, iFq(— p +V; — l) is - strictly speaking 
- undefined and a divergent series. However, (1 —z)"~ p remains well defined as z — ► — 1. 
Therefore, (14.111 ) essentially corresponds to an analytic continuation that implicitly uses the 
concepts of Abel summation which is for instance discussed in Hardy's classic book [78 1. 
Inserting ( 14.111 ) into ( 14. 10b yields: 

G p a \z) = 2P r( P + ° + 1) 1 F 1 (-p;«+l; Z /2). (4.12) 

A confluent hypergeometric series iF\ (a;b;z) converges absolutely for all complex a, b, and 
Z as long as — b ^ No (see for example 11 191 p. 2]). Thus, for — b ^ No such a \F{ is an entire 
function in z. Since we always assume a > — 1 , G p (z) is in every neighborhood of the origin 
z = an analytic function. Thus, ( 14.12b shows that the remarks by Gottlieb and Orszag |49 
p. 42] and by Doha [41 , p. 5452], who stated that such a Laguerre series converges faster than 
algebraically if the function under consideration is analytic at the origin, is imprecise. 

If the argument z of a confluent hypergeometric function i^i (a;b,z) is real and approaches 
+oo, then we have the following asymptotic behavior (see for example 11 1 191 Eq. (4.1.7)]): 

Ma-Mz) - ^jeV- fe [l+0(l/z)], z^+oo. (4.13) 
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This asymptotic estimate shows that the integral 

/V-'z«[G< a) (z)] 2 dz 

= 2P r(Sfir ] f ^^[^o(-P;«+i; Z /2)] 2 ck (4.14) 

converges if a > — 1 and a + 2p > — 1 hold. At the lower integration limit z = 0, the integrand 
behaves like z a , which requires a > — 1 , and at the upper integration limit z = °°, the integrand 
behaves like z~ a ~ 2p ~ 2 , which requires a + 2p > —1. Thus, Gp"'(z) belongs for a > — 1 and 
a + 2p > -1 to the Hilbert space L 2 ae _ z ([0,°°)) defined by (t2TTb. 

It is possible to check the correctness of the power series representation ( 14.12b for Gp"' (z) 
by expanding it in terms of generalized Laguerre polynomials. According to dl.lbb . we then 
have to compute the following inner product: 

jXa) =2 p n\ Y{p + a+\) 



r(a + n + l) r(a + l) 
x / e- z z a Ll a \z) l F 1 (-p;a+l;z/2)dz. (4.15) 



Since the confluent hypergeometric function \F\ in (14.15b is analytic in every neighborhood 
of the origin, integration and summation can be interchanged. We now use [50 Eq. (7.414. 11) 
on p. 850] 

'e-V--e' W d, = gffli±^ , Re, r) >0. (4,6, 



to obtain: 



jT e- z z a Ll a) ( z ) iFj (-p; a + 1; z/2) dz 



T(a + 1) ~ (-p) w (-w)n 2 - m (41?) 



The fact that the Pochhammer symbol (—m)„ with m,n £ No satisfies (— m) n = for m = 
0, 1, ...,« — 1 suggest the substitution m — > « + v, yielding: 



e - z z a 4 a) (z)iFi(-p;a + l;z/2)dz 

= (-')- r(g+ 'l'- p) - ^(-p + -a/2). 



(4.18) 



2"n! 

With the help of dl8l > we obtain 

^(-p+n-,1/2) = (l/2)P-" = 2"-P. (4.19) 
Combination of ( 14.18b and ( 14.19b then yields: 

' e- z z a d a \z) iFj (-p; a + l;z/2) dz 



(-1)" r(a + l)(-p)„ 
2P n! 



(4.20) 
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We also obtain this result if we combine (14.8b and ( 14.15l >. 

As a more complicated example, let us now consider the following Laguerre series: 

2 ,4 a \aMc;z) = t d a \z). (4.21) 

„=0 {C) n (U -f- l)n 

Concerning the real parameters a, b, and c, we assume for the moment only that this series 
exists, which requires — c ^ No, that it does not terminate, which requires —a, — b ^ No. and 
that it converges with respect to the norm ( 12.61 ) of the Hilbert space l^ z a e - z ([0, 00 ))- It follows 
from (12.16b that this is the case if c - a - b + (a+ l)/2 > 0. 

For an investigation of the analyticity of (^^b;c;z) at the origin, let us define 

X ia) = )n (a)n(b)n ^ 

(c)„(a + l)„ 



By proceeding as in the case of the Laguerre series ( 14.7b for G p (z) we obtain after a short 
calculation: 

2 ,s4 a \aMc;z) = £ rt( a + v,b + v;c + v;-l)@^ jj-- . (4.23) 

v=0 {C)vV< l^U + ljv 

It follows from d2~T6l l that T(a + v + n)T(b + v + n)/[T(c + v + «)«!] - n "+b+v-c-\ as „ ^ oo. 
Consequently, the Gaussian hypergeometric series 2 F\ (a + v,b + v;c + v;z) with z = — 1 
diverges for all sufficiently large v € No- Nevertheless, it is possible to associate a finite value 
to the hypergeometric function corresponding to this divergent series with the help of the 
following two linear transformations 111021 p. 47]: 



2Fi(a,b;c;z) = (1 - z)~ a ! 2 *i(a,c -b;c;z/ '(z- 1)) (4.24) 
= (l-z)- b 2Fi{c-a,b;c,z/(z-l)). (4.25) 

The transformation z — > z' = z/(z — 1) maps z = — 1, which is located on the boundary of the 
circle of convergence, to z' = 1/2, which is located in the interior of the circle of convergence. 
Thus, the in general divergent hypergeometric series in ( 14.23b can be replaced by a convergent 
hypergeometric series according to 

2 F l {a + v,b + v-c + v--l) = 2- a - v 2 Fi(a + v,c-b;c + v;l/2) (4.26) 

= 2^ v 2J Fi(c-a,/? + v;c + v;l/2). (4.27) 

Inserting (14.26b and ( 14.27b into (14.23b yields two equivalent power series expansions for the 
function defined by the Laguerre series ( 14.21b . which seem to be new: 

^% Mei z)=2^i^(a + V,c-lr,c + v i l/2)^^^ r (4.28) 

v=0 (C)vV'. {U-\-L) v 



v=0 



A detailed analysis of the domain of analyticity of 2 J^{ (a,b;c;z) requires asymptotic es- 
timates of the behavior of the Gaussian hypergeometric series in ( 14.23b . ( 14.28b and ( 14.291 ) as 
V — > oo. If we use the linear transformation II 1021 p. 47] 

2 F x {a,b;c;z) = (1 -z) c ' a ' b 2 Fi{c-a,c-b;c;z) , (4.30) 
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which certainly holds for \z\ < 1, we obtain: 

2 F l (a + v,b + v,c + v;z) = (1 - z) c ~ a ~ b ~ v 2 F\ (c - a,c - b;c + v;z) . 



(4.31) 



The asymptotic estimate (12.16b yields T(c — a-\-n)T{c— b+n)/[T(c + v+n)n\] ~ yf-a-t>-v-\ 
as n — > oo. Accordingly, the hypergeometric series on the right-hand side of ( 14.311 ) converges 
at least for sufficiently large values of v also for z = — 1 and provides an analytic continuation. 
Moreover, (14.3 lb is a convenient starting point for the construction of an asymptotic expansion 
of the divergent hypergeometric series in ( 14.23b as v — > °° (compare also [ 127, Eq. (15)]): 

2F 1 (a + v,b + v;c + v;-l) = 2^"^ V (-1)" (£Z^Z^ (4 .32) 

„= (c + V)„n! 



2 C 



— a— 6— V 



(c — a)(c — fo) 

(c + v) 



+ 0(v~ 2 ) 



(4.33) 



This asymptotic estimate shows that (a,b;c;z) is for — c £ No analytic in every neigh- 

borhood of the origin. 



5 Exponentially and Factorially Decaying Series Coefficients 

It is immediately obvious that the inner fi series in ( 13.14b converges if the series coefficients 

(cc) 

An decay for sufficiently large indices n exponentially, satisfying for instance 

Ai a) ~ n e R n , eeR, |-R|<1, «^°°. (5.1) 

As the probably most simple example of a Laguerre series of the type of (II. lb with expo- 
nentially decaying coefficients, let us consider the following expansion: 

£ {a) {t;z) = fv4°°(z). (5.2) 

n=0 

It is immediately obvious that the Laguerre series (15.2b converges in the mean with respect to 
the norm (12.6b of the Hilbert space L? ae _-([0,°°)) for |f| < 1, and that it diverges for \t\ > 1. 
Equivalently, we can say that the series ( 12.13b for the norm ||#( a )(?;z)|| a _ z2 produces a 
finite result for \t\ < 1 and that it diverges for \t\ > 1. 

Ignoring for the moment all questions of convergence, let us now assume that t is an 
unspecified complex number and insert = t" into (13.14b . Then, we obtain: 

. i t!f £ (5±l+12t„ (5 . 3) 

v=0 V - n=0 V" 

= E ^^i^o(a + v + l;t). (5.4) 

v=0 V ' 

A nonterminating hypergeometric series \Fo(a;z) converges only in the interior of the unit 
circle. Since we always assume a > — 1, the hypergeometric series iFo^a + V + l;n in ( 15.4b 
does not terminate and we have to require \t | < 1. However, iFq((X + V + l;f) can according 
to ( 13.8b be expressed in closed form, 

^(a + v + ht) = (l-ty a - v -\ reC\{l}, (5.5) 
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which yields an analytic continuation for t j^l. Inserting ( 15.51 ) into d5.4l ) yields 

{-tz/\i-t]y 



g^{t;z) = (1 -t)- a ~ l £ v " /L ; j; (5.6) 

v=0 V - 

= (l-f)- a - 1 e X p(-fz/[l-f]). (5.7) 



Accordingly, the Laguerre series (15.21) for $( a > (t;z) is nothing but the well known generating 
function [TD2] p. 242] 

£4 a) ( Z )^ = (l-rr^rapfc/tf-l]), |f|<l, (5.8) 
«=o 

in disguise. 

As a generalization of the alternating Laguerre series ( 14.71 ). let us consider the following 
expansion with exponentially decaying coefficients: 

(a) _ r(p + q + l) - (-p)„ (g) 

^ p (a,z) - r(a + l) „V S) (a + 1)/" 

peM\N , a + 2p>-\. (5.9) 

Obviously, this Laguerre series converges in the mean for |s| < 1. 

For an investigation of the analyticity of (s;z) at the origin, let us define 

(a) _ r(p + a+l) n-P+n) , y 

A " ~ r(-p) r(a+» + i) ( s) ■ (5 - 1U) 



Inserting this into the modified rearranged Laguerre expansion ( 14.2b yields: 



With the help of ( 13.81 ). we obtain 

iF (-p + v;-s) - (l+*) p ~ v , 5^-1. (5.13) 
Inserting J5.13b into ( 15.121) yields: 

Sf^foz) = (i +j )pE^+1) iFi (_ P; a+ i ; «/[i +J ]). (5.14) 

The argument of the confluent hypergeometric series iFj becomes infinite for s = — 1. If we 
assume a > — 1, then the right-hand side of (15.14b is mathematically meaningful for all z £ C 
and for alls £ C\{ — 1}. However, this is not true for the corresponding Laguerre series ( 15.91 ). 

The convergence of the series ( 12.131 ) for the norm ||Sfp°^(s;z) || a e - z 2 ^ s on ^ guaranteed for 

M < l - 

It follows from (15.141 ) that ^p{s;z) is essentially a confluent hypergeometric series \F\. 
Accordingly, the Laguerre series (E9jfor^ a) ( S ;z) can be rewritten as follows: 

1 F 1 (-p;a+l;«/[l+5]) = (1 + ,)"P f {-sf ±^-L { „ a) (z) . (5.15) 
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If we make in ( 15.151 ) the substitutions p — > — c and s — ► — f, we see that the function Sfp (s;z) 
is nothing but the well known generating function 11091 Eq. (3) on p. 202]) 

t { 1tl\[ X) t " = (l-0- c i^i(^a + l;-*/[l-']) (5-16) 

in disguise. 

If we set s = — 1 in the Laguerre series (15.9b for @p(s;z), we obtain the Laguerre series 

( 13.31 ) for z p . To analyze the behavior of £fp"' (s;z) as s J, —1, we write s = — 1 + 8 with 5 > 
and consider in (15.141 ) the limit 5 j 0. We obtain 

^ o) (*-l;z) = 50 T ( P + ° + 1 } lfl (-p; g + 1; (g - l) z /g) . (5.17) 

Next, we use the following asymptotic estimate Ml 191 Eq. (4.1.8)]: 

iF^a-Mz) = ^(-^[l+oa/z)], z^-oo. (5.18) 

If we insert ( 15.18b into ( 15.171 ). we obtain in the limit of vanishing 5 > 0: 

Sf p (a) (-l;z) = limSf p (a) (5-l;z) 

<5|0 

= ta[(l- 5 ) z ]p{ 1+ o(^)}=A (5.19) 

Thus, ^p w (*;z) 

possesses the one-sided limit z p as s j — 1, as it should according to the 
Laguerre series d3~3l and 

The transformation formula (13.14b can only be used if the function represented by the 
Laguerre series is analytic at the origin, because otherwise the inner jj. series diverges. Thus, 
it is not possible to obtain in this way an explicit expression for a nonanalytic function defined 
by a Laguerre series of the type of ( II. lb with monotone and algebraically decaying series 
coefficients Ai" . 

However, the fact that the in general nonanalytic power function z p can according to (15.19b 
be obtained by considering the one-sided limit s J, — 1 in the analytic function Sfp (s;z) indi- 
cates that the situation is not as hopeless as it may look at first sight. 

(cc) 

Let us therefore assume for the moment that the coefficients X,) of a Laguerre series are 
monotone and decay algebraically as n — > oo. This implies that the inner ji series in (13.14b 
produce infinities. However, we could try to apply the transformation formula ( 13.14b to a 
modified Laguerre series with coefficients t n M . Since these coefficients decay exponen- 
tially for \t\ < 1, this modified Laguerre series represents an analytic function. If we succeed 
in constructing an explicit expression for this analytic function, then we can try to perform the 
one-sided limit 1 1 1 in this expression. If this can be done, we obtain an explicit expression 
for a nonanalytic function. Obviously, this idea deserves to be investigated further. 

Although the explicit expression ( 15.14b for &p a \s;z) is mathematically well defined for 
allseC\{ — 1} and possesses a one-sided limit for s J, — 1 , its Laguerre series (15.9b requires 
more restrictive conditions. In general, generating functions of the kind of ( 15.9b converge 
only for \s\ < 1, i.e., in the interior of the unit circle, and possibly also for some points on the 
boundary of the unit circle. 
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If we make in ( 15.9b or in ( 15.15b the substitution s — ► 1, we obtain the series expansions 
( 14.7b and ( 14.12b . respectively, for Gp"'(x). Thus, in this case it is possible to extend the 
Laguerre series ( 15.91 , which converges exponentially in the mean for \s\ < 1, to the boundary 
of its circle of convergence, yielding a convergent Laguerre series with algebraically decaying 
coefficients that ultimately have strictly alternating signs. 

It is important to note that this approach does not always work. If we make in the La- 
guerre series (15.2b . which because of (15.7b is nothing but the generating function (15.8b . the 
substitution t — » — 1, we formally obtain: 

e z/2 = 2 «+i £ (__i)« L ( a )( z ). (5.20) 

ii=0 

However, e Z//2 does not belong to the Hilbert space i^a e - z ([0, 00 )) defined by ( 12.7b . and it 
follows from (12.18b that the series d2.12b for the norm ||e z / 2 || „ _ z 2 ^ oes not conver g e - 

It is also not possible to assume \s\ > 1 in the Laguerre series <E3 forSf p (a) (>;z), be- 
cause then the series (12.12b for the norm ||^p a ' (s;z) \\ z a e - z 2 diverges. In contrast, the explicit 

expression ( 15.14b for ^p a \s;z) remains meaningful for \s\ > 1. 

If the series coefficients of a power series decay factorially, we can expect that the function 
represented by this power series is analytic in the whole complex plane C. The most obvious 
examples are the exponential exp(z), or - slightly more general - the confluent hypergeo- 
metric series \Fi(a;b;z). In contrast, power series, whose coefficients decay algebraically, 

converge only in compact subset of the complex plane. It should be of interest to study how 

let) 

factorially decaying series coefficients X„ influence the analyticity properties of functions 
represented by a Laguerre series of the type of ( II. lb . 

As an example, let us consider the following Laguerre series: 

«=o \ w + L )n 

Thus, we choose A,!"' = s"/(a + 1)„ in ( II. lb . Inserting this into the modified rearranged 
Laguerre expansion ( 14.2b yields: 

^ a \s;z) = f . ( ~^, )V , T ^ (5.22) 
--^ ( - SZ)V (5.23) 



v =o (a + l) v v! 



The infinite series in ( 15.22b can be expressed as a generalized hypergeometric series 0^1 ■ We 
then obtain: 

^ a \s; Z ) = e s qFi (a + !;-«). (5.24) 



However, this is nothing but a known generating function II 1021 p. 242]. 

The generalized hypergeometric series 0^1 (oc + 1; — sz) in ( 15.24b obviously converges for 
all s,z £ C as long as a + 1 is not a negative integer. Since we always assume a > — 1, we thus 
can conclude that j^ a '(s;z) is an analytic function in the sense of complex analysis in both 
s and z- Moreover, the series (EE) for the norm \\^ a \s;z)\\ z a e - z j converges for arbitrary 
seC. 
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6 Computational Approaches 

The examples considered in Sections H] and [5] show that it is indeed possible to transform a 
Laguerre series of the type of dl.lt with the help of (13 . 14b to a power series about z = of the 
type of (II. 3t . The essential requirement is that the inner fi series in ( 13.141 ) has to converge. 

This is the case if the coefficients A„ (a) of the Laguerre series decay sufficiently rapidly, or 

Ice) 

- if the coefficients A„ only decay algebraically in magnitude - ultimately have strictly 
alternating signs. 

In Section |5] the transformation of Laguerre series with exponentially and factorially de- 
caying series coefficients was discussed. Apart from finite Laguerre expansions, this is pretty 
much the best scenario which can occur in this context. The convergence of the inner fi series 
is guaranteed, and often such a jj. series converges rapidly, making the resulting power series 
expansion computationally useful even if it is not possible to find a convenient closed form 
expression for the inner jj. series. 

As discussed in Section[4] the most challenging and therefore also most interesting prob- 
lems occur if the coefficients A„ decay algebraically as n — ► °°. If algebraically decaying 
coefficients A„ ultimately have the same sign, the inner jj, series in (I3.14t diverges, and it is 
also not possible to sum these series to something finite. Accordingly, the transformation for- 
mula (13 . 14b leads to a power series expansion having infinitely many series coefficients that 
are infinite in magnitude. This simply means that a power series in z does not exist because 
the function under consideration is not analytic at the origin. 

The inner fi series in (13. 14b also does not converge if the coefficients kj,^ decay alge- 
braically and ultimately have strictly alternating signs. But in this case, suitable summation 
techniques are capable of associating something finite to the divergent inner jj. series. 

In all examples considered in Section|4] it was possible to sum the divergent inner fi series 
with ultimately strictly alternating signs by means of explicit analytic continuation formulas 
for hypergeometric series jFq and 2^1, respectively. For example, the Laguerre series ( 14.7b 
for Gp™' (z) leads according to d4.10t to an inner jj, series that can be expressed as a binomial 
series iFq(— p + v; — 1 ) that diverges for fixed p £ R \ No at least for sufficiently large values 
of V € No. Nevertheless, the summation of this divergent series is trivial because every hyper- 
geometric series \Fo(a;z) with z ^ 1 possesses according to (13.8b a very simple closed form 
expression which accomplishes the necessary analytic continuation. 

Similarly, the Laguerre series (14.21b for 2 J%f CC \a,b;c;z) leads according to (14.23b to an 
inner ji series that can be expressed as a hypergeometric series 2^1 {a + V,b + v;c + v; — 1) 
that also diverges for sufficiently large v e No- But again, it is almost trivially simple to 
find explicit analytic continuation formulas that replace this divergent series by convergent 
expansions. The linear transformations ( 14.24b and ( 14.25b do the job. 

Unfortunately, the situation is not nearly as nice if the series coefficients A,, have a more 
complicated structure. Let us for example consider the following Laguerre series with coeffi- 
cients A,, that are ratios of p + 1 numerator and p + 1 denominator Pochhammer symbols: 



Obviously, this Laguerre series generalizes the Laguerre series 

(E2B for 2 -J^i {a, b;c;z). 




peN . 



(6.1) 
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It is easy to show that the application of the transformation formula ( 13.14b to the Laguerre 
series ( 16.1b leads to inner ji series that can be expressed as generalized hypergeometric series 

f at + v,...,a p+1 + v \ _ ^ (ai + v)„---(a p+ i + v)„ z" 
p+ltp \ bl + y,...,b p + v ' Z J ~ L (b l +v) n ---(b p + v) n n\ 

with argument z = — 1. Just like the Gaussian hypergeometric series 2F1, a nonterminat- 
ing generalized hypergeometric series p +\F p converges only in the interior of the unit circle. 
Moreover, ( 12.161 1 implies 

r(a 1+ v + n)r( ap+1 + v + «) ^ n g l+ ... +ap+l+v - bl+ ...- bp n ^ oo 
T{b x +v + n)T{b p + v + n)n\ 

This asymptotic estimate shows that the generalized hypergeometric series ( 16.21 ) with argu- 
ment z = — 1 diverge for sufficiently large values of v and have to be summed. Unfortunately, 
the theory of the generalized hypergeometric series p +\F p with p > 2 is not nearly as highly 
developed as the theory of the Gaussian hypergeometric series 2F1. Many transformation 
and/or analytic continuation formulas for a generalized hypergeometric series p +\F p with 
p > 2 are either not known or at least much more complicated than the corresponding expres- 
sions for a Gaussian hypergeometric series 2^1 (compare also [25 1). We can safely assume 
that our chances of finding convenient explicit analytic continuation formulas, that can ac- 
complish the summation of the divergent generalized hypergeometric series (16.21 ). are rapidly 
approaching zero with increasing p. 

Thus, the rearrangement of Laguerre series via ( 13.14b can only produce closed form ex- 
pressions if algebraically decaying and ultimately strictly alternating series coefficients A„ 
have an exceptionally simple structure, as for example the ones in the Laguerre series ( 14.7b 
and ( 14.21b . In the case of more complicated coefficients X^ a \ the necessary summations of 
the divergent inner ji series will be too difficult to produce closed form expressions of man- 
ageable complexity. 

These and related considerations may tempt a skeptical reader to argue with some justifi- 
cation that all examples considered in Section|4]are actually fairly simple. Therefore, it is by 
no means obvious whether the approach of this article is capable of producing anything new 
beyond a rederivation of known and comparatively simple formulas. 

While the first conclusion is certainly correct, the second is in my opinion overly pes- 
simistic. The reason is that summation via analytic continuation formulas is not the only 
possibility: We can also try to use numerical summation techniques. Of course, knowing only 
numerical approximations to the leading coefficients of a power series is not nearly as nice 
as knowing explicit and possibly even simple expressions for the series coefficients, but it is 
certainly better than nothing. In practical applications, this limited information may suffice. 

The use of numerical summation techniques is by no means a new idea. Already in 1882, 
Holder ll80l proposed - inspired by an article by Frobenius flTll - to determine the value 
of a power series on the boundary of its circle of convergence with the help of numerical 
summation processes based on weighted arithmetic means. 

As discussed in Appendix [Bj the numerical processes proposed by Holder ll80ll and oth- 
ers subsequently evolved to a sophisticated mathematical theory of so-called regular matrix 
transformations that have many advantageous theoretical features and that can be used for 
the summation of divergent alternating series possessing similar divergence properties as the 
generalized hypergeometric series ( 16.2b . 

Unfortunately, regular matrix transformations are in general at best moderately powerful. 
Therefore, I prefer to use instead nonlinear sequence transformations, whose theoretical prop- 
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erties are not yet completely understood, but which are according to experience often much 
better suited to achieve highly accurate summation results. 

It may be interesting to note that summation techniques such as conformal mappings, re- 
expansions, and Pade approximants were discussed in an article by Skorokhodov [ 1 16 1 on the 
analytic continuation of a divergent generalized hypergeometric series p +\F p with argument 
\z\ > 1. Skorokhodov completely ignored the theoretically much simpler, but also less power- 
ful regular matrix transformations. Nevertheless, I do not think that the summation techniques 
considered by Skorokhodov give best results in the case of the divergent, but summable inner 
ji series occurring in this article. In my opinion, (much) better results can be obtained with 
the help of those nonlinear sequence transformations that are reviewed in AppendixlBl 

In this article, I apply as numerical summation techniques Wynn's epsilon algorithm 
dB.4l . and the two Levin-type transformations ( IB. 14b and ( IB. 151 ). Wynn's epsilon algo- 
rithm produces Pade approximants according to ( IB. 51 ) if the input data are the partial sums 
fn(z) = Li=o7<:^ °f a power series. The two Levin-type transformations ( IB. 14b and ( IB. 15b 
are based on the remainder estimate ( IB. 131 ) which corresponds to the first term neglected in 
the partial sum. Both are known to be highly effective in the case of both convergent and 
divergent alternating series. 

The discussion in Section @] should provide convincing evidence that for z = — 1 and in 
particular for large values of v € No it is not a particularly good idea to use the hypergeometric 
series 2.F1 (a + V , + V ; c + v ; z) for the evaluation of the hypergeometric function it represents . 
With the help of analytic continuation formulas like (14.24b and (14.25b or also (14.30b , com- 
putationally much more convenient hypergeometric series can be derived. However, these 
simplifications are only possible in the case of a Gaussian hypergeometric series 2^1, but 
not necessarily in the case of more complicated generalized hypergeometric series p +\F p , let 
alone in the case of divergent, but summable inner fi series that result from purely numerical 
Laguerre series coefficients A„ . Thus, the hypergeometric series 2^1 (a + V,b + v;c + v; — 1) 
serves as a model problem for other divergent alternating series whose terms also increase in 
magnitude like a fixed power of the index. It should be interesting to see how much can be 
accomplished by employing powerful nonlinear sequence transformations. 

In Tables[T|and|2] the nonlinear transformations mentioned above are applied to the partial 
sums 

s n = s n{ a M = ii-lt^±^^ (6.4) 

of the Gaussian hypergeometric series 2^1 {a + V,b + v;c + v; — 1) with a = 3/2, b = 7/3, and 
c = 21/4. In TableQ]we have v = 0, and in Table|2]we have v = 10. All transformation results 
of this article were obtained using the floating point arithmetics of Maple 11, and the "exact" 
result in Tables [T] and [2] was produced by the Maple procedure hypergeom which computes 
generalized hypergeometric series. 

The entries in the third column of Table Q] are chosen according (IB. 8b . and the entries in 
columns 4 and 5 according to ( IB. 16b . Since Wynn's epsilon algorithm produces according to 
(IB. 5b Pade approximants if the input data are the partial sums of a power series, the entries in 
the third column can be identified with the staircase sequence [L(» + 1)/2J/L«/2J] of Pade 
approximants to the hypergeometric series 2F1 (a + v,b + v;c + v;z) with z = — 1. Here, |xj 
stands for the integral part of x, which is the largest integer v satisfying v < x. 

The asymptotic estimate d2.16b shows that the terms (— l)"(a) n (b) n /[c n nl] of the hyperge- 
ometric series in TableQ]decay in magnitude like n~ 29 / 12 as n — > °°. Therefore, this hyperge- 
ometric series converges, albeit quite slowly. This assessment is confirmed by the data in the 
second column of Table [T] The transformation results in Table [T]indicate that all three trans- 
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Table 1: Application of Wynn's epsilon algorithm and the Levin-type transformations 
d^ifijSn) and 5^ (/?,£„) with j3 = 1 to the partial sums of the slowly convergent hyper- 
geometric series 2 F\ (a + V,b + v;c + v; -1) with a = 3/2, b = 7/3, c = 21/4, and v = 0. 



n 


s„ fz) 


>-2L»/2J) 
2|«/2J 


d 0) (l sa(z)) 


5„ (0) fl sate)) 




Eq. $6A\ 


Eq. dB~4b 


Eq. (|B.14|i 


Eq. (|B.15|l 


u 






i nnnnnnnnnnnnnnn 


1 nnnnnnnnnnnnnnn 


i 


M ^^^^^^^^ 


^^^^^^^^^^^^^^^ 


n f>r\rtr\r\r\r\r\r\r\r\r\r\r\r\ 
u . ouuuuuuuuuuuuuu 


v.ovvvvvvvvvvvvvv 


z 


U. ////// / 


U.uUUUUUUUUUUUUUU 


u.jyou/y/j/io/i'uy 


u. jyou / yZjz 1 s/yoy 


3 


0.46785866 


0.595184349134688 


0.597361776678518 


0.597283629053888 


4 


0.69325438 


0.597114931459132 


0.597128156362257 


0.597151427941970 


5 


0.52349689 


0.597142098567676 


0.597159499152400 


0.597156391214500 


6 


0.65507045 


0.598028964909470 


0.597156082629105 


0.597156376752544 


7 


0.55064699 


0.597156986187236 


0.597156397164508 


0.597156374043337 


8 


0.63518026 


0.597156454880736 


0.597156372412882 


0.597156373980610 


9 


0.56559242 


0.597156266466529 


0.597156374069016 


0.597156373980877 


10 


0.62370437 


0.597156373786525 


0.597156373977124 


0.597156373980968 


11 


0.57457047 


0.597156373530068 


0.597156373981079 


0.597156373980973 


exact 




0.597156373980973 


0.597156373980973 


0.597156373980973 



formations are able to accelerate the convergence of the hypergeometric series 2^1 effectively. 
The least effective, but still very powerful accelerator is Wynn's epsilon algorithm, and the 
most effective transformation is the Levin-type delta transformation which is - as documented 
by numerous references mentioned in Appendix 151- known to be highly effective in the case 
of both convergent and divergent alternating series. 

In Table |2j only Wynn's epsilon algorithm and the delta transformation are displayed. 
The data in the second column show that the partial sums ( 16.41 ) with v = 10 diverge rapidly. 
Nevertheless, it is possible to obtain highly accurate summation results. As in Table Q] the 
delta transformation was clearly more effective than Wynn's epsilon algorithm. 

It is in my opinion remarkable that Levin's transformation ( IB. 14b . whose action on the 
partial sums ( 16.4b is not displayed in Table [2] turned out to be only roughly as effective as 
Wynn's epsilon algorithm. Under the same conditions as in Table|2j I obtained the following 
summation results: 

dfj(l,s (z)) = 0.142894441246992 x 10" 2 , (6.5) 
d^(l,s (z)) = 0.142895355463039 x 10~ 2 . (6.6) 

In view of the in general very good reputation of Levin's sequence transformation, this com- 
paratively weak performance is somewhat puzzling. 

In the case of more complicated generalized hypergeometric series p +iF p , the same gen- 
eral pattern was observed. If the same nonlinear sequence transformations as in Tables [T] and 
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Table 2: Summation of the hypergeometric series 2,F\ (a + v, b + v; c + v; — 1 ) with a = 3/2, 
b = 7/3, c = 21/4, and v = 10 with the help of Wynn's epsilon algorithm and the Levin-type 
transformation (P,s„) with j3 = 1. 



*(*) 5i 0) (l,, (z)) 

Eq. dEl Eq. (El Eq. drTBl) 

0.10000 x 10 1 0.100000000000000 x 10+ 1 0.100000000000000 x 10+ 1 

1 -0.83005 x 10 1 -0.830054644808743 x 10+ 1 -0.517662391110501 x 10+° 

2 0.39395 x 10 2 -0.517662391110501 x 10+° 0.355800262535272 x 10+° 

3 -0.13894 x 10 3 0.176355932965950 x 10+ 1 -0.979095090620413 x 10"' 

4 0.40421 x 10 3 0.150567485119558 x 10+° 0.156275192732603 x 10-' 

5 -0.10245 x 10 4 -0.292708508243714 x 10+° 0.539203120343906 x 10~ 3 

6 0.23385 x 10 4 -0.235266825737199 x 10" 1 0.137225149988457 x 10~ 2 

7 -0.49149 x 10 4 0.351057409892771 x 10-' 0.145137993614918 x 10~ 2 

8 0.96599 x 10 4 0.394766824723114 x 10~ 2 0.142550654190640 x 10~ 2 

9 -0.17957 x 10 5 -0.113107334634727 x 10~ 2 0.142932832720914 x 10~ 2 

10 0.31849 x 10 s 0.127361316374322 x 10~ 2 0.142892054773978 x 10~ 2 

11 -0.54255 x 10 5 0.155543251021275 x 10~ 2 0.142895407283313 x 10~ 2 

12 0.89251 x 10 s 0.143471035609870 x 10~ 2 0.142895196573040 x 10~ 2 

13 -0.14240 x 10 6 0.142504106397527 x 10~ 2 0.142895206303027 x 10~ 2 

14 0.22113 x 10 6 0.142883205557579 x 10~ 2 0.142895205998921 x 10~ 2 

15 -0.33524 x 10 6 0.142902170361815 x 10~ 2 0.142895206004165 x 10~ 2 

16 0.49740 x 10 6 0.142895322013589 x 10~ 2 0.142895206004153 x 10~ 2 

17 -0.72381 x 10 6 0.142895148436317 x 10~ 2 0.142895206004152 x 10~ 2 

exact 0.142895206004152 x 10~ 2 0.142895206004152 x 10~ 2 



|2] are applied to the partial sums 



of the generalized hypergeometric series 3^2(0 + V,b + v,c + v;d + v,e + v; — 1) with a = 3/2, 
b = 7/3, c = 11/5, d = 22/7, and e = 32/1 1, we obtain for the convergent hypergeometric 
series with v = 
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and for the divergent hypergeometric series with v = 10 

= 0.816454762672306 x 1(T 3 , (6.12) 

d$(l,s {z)) = 0.816459448502108 x 10~ 3 , (6.13) 

8$(l,s (z)) = 0.816458731770118 x 10 -3 , (6.14) 

hypergeom = 0.816458731770118 x 10~ 3 . (6.15) 

In the case of the partial sums 

sn = s n (a,b,c,d,e,f,gv) = £ -if '\ '\ ' y (6.16) 

k~0 {e + v)k{f + v) k {g + v) k k\ 

of the generalized hypergeometric series 4^3(0 + V,b + v,c + v,d + v;e + v,f + V,g + v; — 1) 
with a = 3/2, b = 7/3, c = 11/5, d = 16/17, e = 18/19, / = 22/7, and g = 32/1 1, we obtain 
for the convergent hypergeometric series with v = 

e$ = 0.538509188429164, (6.17) 

<*i ( 2 (Mo(z)) = 0.538509188330837, (6.18) 

8$(l,s (z)) = 0.538509188330835, (6.19) 
hypergeom = 0.538509188330835, (6.20) 

and for the divergent hypergeometric series with v = 10 

e$ = 0.819691553980977 x 10" 3 , (6.21) 

d$(l,so(z)) = 0.819696193819277 x 10 -3 , (6.22) 

8$(l,s (z)) = 0.819695479036364 x 10 -3 , (6.23) 

hypergeom = 0.819695479036364 x 10" 3 . (6.24) 

It was emphasized in Section|5]that the inner fi series in ( 13.141 i converge if the coefficients 

(cc) 

An of a Laguerre series decay exponentially as n — > 00. Accordingly, summation techniques 
are not needed in this case. Nevertheless, nonlinear sequence transformations can be ex- 
tremely useful even in the case of exponentially decaying Laguerre series coefficients. Let us 
for example assume that the coefficients A,!"' of a Laguerre series satisfy 

= t"il" ] (6.25) 

and that it is not possible to derive a closed form expression for the inner ji series in ( 13.141 ). 
Thus, the the inner pi series have to be evaluated numerically. 

(n) 

Let us now also assume that the coefficients t„ increase like a fixed power of n as n — > 00. 
In spite of this unfavorable behavior, the inner ji series converge as long as \t\ < 1, However, 
convergence can become prohibitively slow if |f | is only slightly smaller than one. The con- 
vergence problems are particularly severe if the terms in inner pi series ultimately have the 
same sign because then the fi series do not converge for t = 1 and are also not summable. 
But again, nonlinear sequence transformations can be extremely useful to speed up the con- 
vergence of such a monotone series with (very) slowly decaying terms (see for example [85| 
or ll27l Section 2.2.6] and references therein). 
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As documented by the recent books by Cuyt, Brevik Petersen, Verdonk, Waadeland, and 
Jones [36] and by Gil, Segura, and Temme or by a review by Temme [128|, there is 
currently a lot of work being done on the efficient and reliable evaluation of special functions. 
As I had shown in several articles iU [85] [J36j [J37] [E3 [140] [143] Mj\ [H3 ESI [1621 . 
nonlinear sequence transformations can be extremely useful in this respect. 



7 Guseinov's Rearranged One-Range Addition Theorems 

The analyticity of Laguerre series is a problem of classical analysis, but I became interested 
in this mathematical topic because of some open questions in molecular electronic structure 
theory. During the work for my forthcoming review on addition theorems [ 1 54 1 , 1 came across 
some articles by Guseinov ll53l l54l l56l who had constructed one-range addition theorems for 
Slater-type functions 1117111181 . In unnormalized form, Slater-type functions, which play a 
major role as basis functions in atomic and molecular electronic structure calculations, are 
defined as follows: 

Z ^(j3,r) = (Pr) N - L - l e-P^ L M (lir)- (7-D 

Here, r G R 3 , ^ M (/3r) = {fir) L Y^ (d ,<$>) is a regular solid harmonic and T^(0,0) is a (sur- 
face) spherical harmonic, j3 > is a scaling parameter, N G R \ N is a kind of generalized 
principal quantum number which is often, but not always a positive integer > L + 1 , and L and 
M are the usual (orbital) angular momentum quantum numbers. 

Let us assume that {<p^(r)}„ ^ m is a complete and orthonormal function set in the Hilbert 
space 



L 2 ( 



5 ) = {/:R 3 ^c|| |/(r)| 



|2 d 3 r<ocl (7.2) 



of functions that are square integrable with respect to an integration over the whole R 3 . 
Since any / G L 2 (R 3 ) can be expanded in terms of the complete and orthonormal functions 
{ ( P™e( r )}n,e,m, a one-range addition theorem for f(r ± r') can be formulated as follows: 

f(r ± r') = £ ±r') < f (r) , (7.3a) 

nim 

C:/f;±r>) = J [<,(r)]7(r±r')d 3 r. (7.3b) 

The expansion ( 17.31 ), which converges in the mean with respect to the norm of the Hilbert 
space L 2 (R 3 ), is a one-range addition theorem since the variables r and r' are completely 
separated: The dependence on r is entirely contained in the functions <p™» (r), whereas r' 
occurs only in the expansion coefficients C™»(/;±r') which are overlap or convolution-type 
integrals. 

One-range addition theorems of the kind of (17.3b were constructed by Filter and Steinborn 
B31 Eqs. (5.11) and (5.12)] and later applied by Kranz and Steinborn |98| and by Trivedi 
and Steinborn [132]. An alternative derivation of these addition theorems based on Fourier 
transformation combined with weakly convergent expansions of the plane wave exp(±ip • r) 
with p,r G R 3 was presented in 11351 Section VII]. 

As discussed in 11511 Section 3]), it is also possible to formulate one-range addition theo- 
rems that converge with respect to the norm of a weighted Hilbert space 

L 2 (R 3 ) = {/: R 3 -> C || w(r) |/(r)| 2 d 3 r < «.} , (7.4) 
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where w(r) ^ 1 is a suitable positive weight function. If we assume that / e Ll(R 3 ) and 
that the functions { VC/( r )}n tjn are complete and orthonormal in L^,(R 3 ), then we obtain the 
following one-range addition theorem 01511 Eq. (3.6)]): 



f(r ± r') = £ CUf, w; ±r') <>(r) , (7.5a) 



C^(/,w;±r') = / [<>(r)]*w(r)/(r±r')d 3 r. (7.5b) 



A one-range addition theorem for a function / : R 3 — > C is a mapping R 3 x R 3 — > C. Com- 
pared to the better known two-range addition theorems like the so-called Laplace expansion 
of the Coulomb or Newton potential 1/r, which possesses a characteristic two-range form 
(see for example 01511 Eq. (1.2)]), one-range addition theorems have the highly advantageous 
feature that they provide a unique infinite series representation of f(r ± r') with separated 
variables r and r' that is valid for the whole argument set R 3 x R 3 . Further properties of ad- 
dition theorems in general and of one-range addition theorems in particular will be discussed 
in my forthcoming review [154]. 

In his one-range addition theorems, Guseinov used as a complete and orthonormal func- 
tion set the following functions 1531 Eq. (1)], which - if the mathematical notation for the 
generalized Laguerre polynomials is used - can be expressed as follows B1511 Eq. (4.16)]: 



(n + £ + k+l)! J n - 1 - 1 V ' ' e V ' 

neN, k= -1,0,1,2,... , 7>0. (7.6) 



As discussed in the text following [ 151 , Eq. (4.20)], Guseinov's functions d7.61 > can - depend- 
ing on the value of A: = — 1,0, 1,2,... - reproduce several other physically relevant complete 
and orthonormal function sets. 

Guseinov's functions are orthonormal with respect to the weight function w(r ) = r k (com- 
pare also EH Eq. (4)]): 

/ K*(y, r )] * rk iX/Ar* r ) dV = 5 ""' s w • ( ? - ? ) 

Accordingly, Guseinov's functions are complete and orthonormal in the weighted Hilbert 
space 



L*(K 3 ) = {/:R 3 ^C|| r*|/(r)| 2 d 3 



V<oo| ) *:= -1,1,2,... . (7.8) 

For k = 0, we retrieve the Hilbert space L 2 (R 3 ) of square integrable functions defined by ( 17.21 ). 

As long as the principal quantum number N is not too negative, a Slater-type function 
Xnl(P i r ) i s f° r a fixed value of k an element of the weighted Hilbert space L 2 ^ (R 3 ). In this 
case, Guseinov's approach J53j|54][56), who constructed one-range addition theorems by ex- 
panding %^ L {^,r ±r') in terms of his complete and orthonormal functions { k ^"^(y, r )} n ,f,m 
with in general different scaling parameters j3 ^ y > 0, is mathematically sound. For fixed 
£=—1,0,1,2,..., Guseinov constructed expansions that converge in the mean with respect 
to the normofL^(R 3 ). 

However, Guseinov replaced in his one-range addition theorems his complete and or- 
thonormal functions k X I"?t(P,i') by nonorthogonal Slater-type functions with integral princi- 
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pal quantum numbers via 11511 Eq. (6.4)] 



(2j3) k + 3 (n + £ + k+ 1)P 1/2 



(n-t-l)l 

x £fa (2e+k+v+2)\ v \ Xv+MA ^ r >- (7 ' 9) 

This is still legitimate. However, Guseinov also rearranged the order of summations of the 
resulting expansions. In this way, Guseinov formally constructed expansions of Slater-type 
functions j^f L (j3,r ±r') with in general nonintegral principal quantum numbers N e R\N 
in terms of Slater-type functions XnAP' r ) w i m integral principal quantum numbers n 6 N 
located at a different center (see also 111511 Section 6]). 

Slater-type functions are complete in all Hilbert spaces, which Guseinov implicitly used 
(for an explicit proof, see |9"T1 Section 4]), but not orthogonal. Thus, Guseinov's approach 
corresponds to the transformation of an expansion in terms of a complete and orthogonal 
function set to an expansion in terms a of a complete, but nonorthogonal function set. 

Unfortunately, the completeness of a nonorthogonal function set in a Hilbert space does 
not suffice to guarantee that an essentially arbitrary element of this Hilbert space can be ex- 
panded in terms of this function set (the nonanalyticity of certain Laguerre series discussed 
in this article is just another confirmation of a much more general fact). This insufficiency 
is well documented both in the mathematical literature (see for example 11381 Theorem 10 
on p. 54] or lf79l Section 1.4]) as well as in the literature on electronic structure calculations 
ll88l [89] |90] HU [92] [93)), but nevertheless frequently overlooked. Horrifying examples of 
nonorthogonal expansions with pathological properties can be found in |[89l Section III.I]. 

Consequently, it is not at all obvious whether the mathematical manipulations, that pro- 
duced Guseinov's rearranged one-range addition theorems, are legitimate and lead to expan- 
sions that are mathematically meaningful. This has to be checked. So far, Guseinov has cat- 
egorically denied that there might be any problem with the legitimacy of his rearrangements 
ED PP- 8 and 23 -24]. 

Since Guseinov's original addition theorems are expansions in terms of generalized La- 
guerre polynomials, and since the radial parts of Slater-type functions are after the cancella- 
tion of common exponentials nothing but powers, the results of this article about the analytic - 
ity of Laguerre series could in principle be used to investigate whether Guseinov's rearrange- 
ments are legitimate. For that purpose, it would be necessary to determine the decay rates 
and the sign patterns of the coefficients of the generalized Laguerre polynomials occurring in 
Guseinov's original addition theorems. 

Unfortunately, one-range addition theorems for Slater-type functions are fairly compli- 
cated mathematical objects, and the coefficients of the generalized Laguerre polynomials are 
according to ( 17.5b essentially three-dimensional overlap integrals. Thus, we would be con- 
fronted with enormous and possibly even unsurmountable technical problems if we try to 
analyze the decay rates and sign patters of the coefficients of these Laguerre series. 

Fortunately, some insight can be gained by analyzing not the comparatively complicated 
one-range addition theorems, but their much simpler one-center limits. These are also expan- 
sions in terms of generalized Laguerre polynomials, albeit with much simpler coefficients. It 
will become clear later that this approach cannot answer all questions in interest. Neverthe- 
less, it is better than nothing since we can obtain at least some nontrivial answers. 

Let us now consider the one-center expansion of a Slater-type function Xnl (fi > r ) w i tri m 
general nonintegral principal quantum numbers N G R \ N in terms of Guseinov's functions 
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v ^'n{(P' r ) w i*h equal scaling parameters j3 > II 1521 Eq. (5.7)]: 
X$t(P,r) = ( V 2N _ X T(N + L + k + 2) 



x 2 



(~N+L+l) v M ,o. 



1/2 



o [(v + 2L + fc + 2)!v!] 
JVeR\N, j3>0, fc = -1,0,1,2,... . (7.10) 

If 7Y G N and 7Y > L + 1, the infinite series on the right-hand side terminates because of the 
Pochhammer symbol ( — N + L + 1 ) v . 

The expansion ( 17.101 ). which corresponds to the one-center limit r' = of Guseinov's 
one-range addition theorem for ^^f L (j3,r ± r') with equal scaling parameters, is a special 
case of the Laguerre series (13.31 ) for z p with in general nonintegral p G M\No- As discussed 
in Section [3] a rearrangement of the Laguerre series (13.31) for z p is legitimate if and only if 
p is a nonnegative integer, i.e., if p = m € No- If p ^ No, z p is not analytic at z = and we 
formally obtain the power series ( 13.71 ). which is not a mathematically meaningful object since 
an infinite number of its power series coefficients are according to d3.9l l infinite in magnitude. 

Thus, in the case of equal scaling parameters /3 > 0, the one-center limit r' = of Gu- 
seinov's rearranged addition theorem for Slater-type functions Xnl(P ,i"±r') does not exist 
if the principal quantum number is nonintegral, i.e., if N e M \ N. 

Let us now consider the one-center expansion of a Slater-type function Xnl(P i r ) w i tn m 
general nonintegral principal quantum numbers N 6 R \ N in terms of Guseinov's functions 
k^n e(Y> r ) witri different scaling parameters j3 ^ y > 11521 Eq. (5.9)]: 



(2y)LHk+3)/2pN-l r(N + L + k + 2) 



[p + Y }N+L+k+2 (2L + k + 2 
V + 2L + k + 2) 



xE 

v=0 



V 



1/2 

k^V+L+l,L(Y: r ) 



x 2 Fi ^-v,N + L + k + 2-2L + k + 3;^-^j . (7.11) 

The expansion ( 17.111 ). which corresponds to the one-center limit r' = of Guseinov's one- 
range addition theorem for Q3 , r ± r 1 ) with different scaling parameters, is a special case 
of the following Laguerre series H151I Eq. (6.12)]: 

zPe- = (i- M r-p- ir(a+p+1) 



na + l) 



X 

n=0 



£ 2^1 ^-n,a + p + l;a+l;^j 4 a) (z), 
peM\N , Re(p + a)>-l, ue(-»,l/2). (7.12) 



This expansion can be derived with the help of ( 13.21 ). The condition — °° < u < 1/2 is necessary 
to guarantee its convergence in the mean with respect to the norm ( 12.6b of the weighted Hilbert 
space L 2 „ e _-([0,oo)). For u = 0, the terminating Gaussian hypergeometric series 2^1 can be 
expressed in closed form with the help of Gauss' summation theorem II 1021 p. 40], and ( 17.121 ) 
simplifies to give (13.3b . 

If we insert the explicit expression ( 12.11 ) of the generalized Laguerre polynomials into 
(17 . 1 2b and interchange the order of summations, we also obtain a formal power series in z. 
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Unfortunately, an analysis of the resulting power series becomes very difficult because of the 
terminating Gaussian hypergeometric series %F\ in ( 17.121 1. An analysis of the behavior of this 
2F1 asn^M would most likely be a nontrivial research project in its own right. However, 
we can argue that z p exp(uz) is only analytic at z = if p = m € No, yielding the expansion 
z"'exp(uz) — L"=o u n z m+n /n\. If p is nonintegral, a power series expansion of z p exp(uz) 
about z = does not exist. 

Thus, also for different scaling parameters j3 ^ 7, the one-center limit r' = of the rear- 
ranged addition theorems for Xnl(P ■,i'±r / ) does not exist if the principal quantum number 
N is nonintegral, i.e., if JV e M \ N. 

These observations are quite consequential: Rearranged one-range addition theorems for 
Slater-type functions Xjfr{fi,r±r') with nonintegral principal quantum numbers N 7^ 1,2,... 
play a central role in numerous articles by Guseinov and coworkers on one-range addition 

theorems |^|^^^H^^|^|^EE|^|^^CD1I23EZ3[S1C3I23[Z2I. 

By analyzing the rearrangement of the Laguerre series of the comparatively simple func- 
tions z p and z p ex.p(uz), we could arrive at some conclusions about the legitimacy of Gu- 
seinov's approach. Nevertheless, some interesting questions are still open. For example, 
the nonanalyticity arguments presented here allow no conclusions about the validity of Gu- 
seinov's rearrangements in the case of Slater-type functions with integral principal quantum 
numbers. 

Another interesting but open question is whether Guseinov's rearrangements produce in 
the case of nonintegral principal quantum numbers one-range addition theorems that are in- 
valid for the whole argument set R 3 x R 3 , or whether only the one-center limits of these 
addition theorems are invalid. This is a practically very relevant question. If only the one- 
center limits are invalid, then it would be conceivable that Guseinov's rearranged one-range 
addition theorems might be mathematically meaningful or possibly even numerically useful 
in a restricted sense as approximations, although they do not exist for the whole argument 
set R 3 x R 3 . This remains to be investigated. But the burden of proof lies in all cases with 
Guseinov. 

8 Summary and Conclusions 

The generalized Laguerre polynomials belong to the so-called classical orthogonal polynomi- 
als of mathematical physics, and they are characterized by the orthogonality relationship ( 12.4b 
that involves an integration over the semi-infinite positive real axis. Accordingly, generalized 
Laguerre polynomials can be used for the representation of functions on unbounded domains, 
and in particular also for the representation of the radial parts of functions /: R 3 — > C ex- 
pressed in terms of spherical polar coordinates. 

It is generally accepted that expansions in terms of orthogonal polynomials have many 
highly advantageous features. However, there is one undeniable drawback: Normally, ex- 
pansions in terms of orthogonal polynomials converge in the mean with respect to the norm 
of the corresponding Hilbert space, but not necessarily pointwise. Accordingly, orthogonal 
expansions are not necessarily the best choice if the local properties of a function matter. 
As documented by the popularity of Pade approximants, power series are also not free of 
weaknesses, but at least in the vicinity of the expansion point, power series are normally very 
convenient and very useful for an accurate description of the local properties of a function. 
Therefore, the construction of power series expansions for functions defined by orthogonal 
expansions should be of principal interest. 

This article describes a complementary treatment of Laguerre series of the type of ( II. 11 1. 
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Normally, one starts from a known function f(z) belonging to the weighted Hilbert space 

L^a e - : ([O7 00 )) defined by ( 12.7b . and one tries to determine the coefficients A« via ( 1 1 . lbb by 
exploiting the orthogonality of the generalized Laguerre polynomials. 

(cc) 

In this article, it is instead assumed that only the Laguerre series coefficients X,) are 
known, either in the form of explicit expressions or numerically, but not the function /(z). 
With the help of the transformation formula ( 13.141 ). it is then possible to construct a formal 
power series expansion of the unknown function represented by the Laguerre series. 

This approach does not guarantee success since there are many functions which belong 
to the Hilbert space L^, e _ z ([0,°°)) but which are not analytic at the origin. A simple ex- 
ample of such a nonanalytic function possessing a Laguerre series is the power function z p 
with nonintegral p € R\No- As discussed in Section [5] it is possible to construct a formal 
power series for z p from its Laguerre series (13.3l l. but the resulting power series is mathemat- 
ically meaningless since it contains infinitely many power series coefficients that are infinite 
in magnitude. 

Thus, the key question is whether the power series for the unknown function obtained 
via ( 13.141 ) is mathematically meaningless, or whether this power series represents an analytic 
function in the sense of complex analysis. It seems that this question has not been treated 
properly in the literature yet. I am only aware of short remarks by Gottlieb and Orszag |49, 
p. 42] and by Doha ATI p. 5452], respectively, who had stated that a Laguerre series of the 
type of (11.11 ) converges faster than algebraically if the function under consideration is analytic 
at the origin. But this statement is imprecise and ignores the pivotal role played by divergent, 
but summable inner ji series in ( 13.141 ). Some general aspects of the summation of divergent 
series are reviewed in AppendixlAl 

By analyzing the convergence properties of the inner jj. series in (13.141 ). some simple 
sufficient conditions can be formulated which guarantee that the resulting power series is 
mathematically meaningful and represents an analytic function. 

As discussed in Section [5] the most benign situation occurs if the Laguerre series coef- 

(oc) 

ficients A„ decay exponentially or even factorially as n — > 00. Then, the inner fi series in 
( 13.14b converge and the resulting power series is mathematically meaningful and represents 
an analytic function. In this way, numerous generating functions for the generalized Laguerre 
polynomials can be rederived easily. 

As discussed in Section 01 a much more interesting situation occurs if the coefficients 

(cc) fee) 

An decay algebraically in magnitude as n — > 00. If algebraically decaying coefficients A„ 
ultimately have the same sign, the inner ji series in ( 13.141 ) diverge. This alone wound not nec- 
essarily be such a bad thing, but it is not possible to sum these divergent series to something 
finite. Accordingly, the transformation formula ( 13.14b leads to a power series expansion hav- 
ing infinitely many series coefficients that are infinite in magnitude. This simply means that a 
power series of the type of (11.3b does not exist because the function under consideration is not 
analytic at the origin. An example is the Laguerre series (13.3b for z p : Its series coefficients 
decay algebraically in magnitude and ultimately all have the same sign. This observation 
suffices to show once more that z p is not analytic at the origin if p 6 R \ No- 

(cc) 

The situation changes radically if the coefficients A„ decay algebraically in magnitude 
as n — > °°, but ultimately have strictly alternating signs. Then, the inner fi series in (13.14b 
still do not converge, but now summation techniques can be used to associate finite values to 
these divergent alternating series. In such a case, the formal power series obtained via ( 13.14b 
is mathematically meaningful and represents an analytic function. 

For example, the additional sign factor (—1)" introduced into the Laguerre series ( 13.3b 
for z p yields the Laguerre series ( 14.7b . The application of (13.14b to (14.71 ) leads to divergent 
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inner /I series which correspond to divergent hypergeometric series \Fq. But these series can 
be summed by analytic continuation since they are special cases of the binomial series ( 13.8b . 
Thus, the transformation formula ( 13.14b ultimately produces according to ( 14.121 ) a confluent 
hypergeometric series \F\, which is an analytic function in every neighborhood of the origin 
and which is also a known generating function of the generalized Laguerre polynomials. 

The summability approach pursued in Section H] can also be used in the case of Laguerre 
series with more complicated coefficients. An example is the Laguerre series (14.211) . The 
application of the transformation formula ( 13.14b leads to divergent, but summable fi series 
that can be expressed as a divergent Gaussian hypergeometric series ^Fi. But again, the 
summation of this divergent series is almost trivially simple because many convenient analytic 
continuation formulas for a ^Fi are known. 

Thus, the summation of a divergent inner ji series can be accomplished by explicit analytic 
continuation formulas if it can be expressed either as a divergent binomial series \Fq or as a 
divergent Gaussian hypergeometric series ^Fi ■ Unfortunately, this is no longer possible if 
the ji series corresponds to a divergent generalized hypergeometric series p +\F p with p > 2. 
As discussed in Section [6] analytic continuation formulas for more complicated generalized 
hypergeometric series are either not known at all or at least much more complicated than the 
corresponding formulas for Gaussian hypergeometric series. 

Thus, we can only hope to find convenient analytical expressions for power series coef- 
ficients Yn if the corresponding algebraically decaying and ultimately strictly alternating La- 

( ex) 

guerre series coefficients A„ possess a very simple structure. But this is a typical limitation 
of all analytical manipulations. As a viable alternative, we can try to use instead techniques 
that accomplish a summation of divergent inner fi series by purely numerical means. Such an 
approach has the additional advantage that it can be applied if only the numerical values of a 
finite set of Laguerre series coefficients xji are available. 

In Section[6] it is shown that certain nonlinear sequence transformations, whose properties 
are reviewed in Appendix [B] are indeed able to sums divergent alternating hypergeometric 
series p +\F p with p > 1 quite effectively. It may be surprising for nonspecialists that the best 
summation results were not obtained by Wynn's celebrated epsilon algorithm (IB. 4b . which 
produces Pade approximants if the input data are the partial sums of a power series, but by the 
Levin-type transformation ( IB. 15b . 

In Section [7] the legitimacy of the rearrangement of Guseinov's one-range addition the- 
orems for Slater-type functions with in general nonintegral principal quantum numbers is 
analyzed. Originally, Guseinov ll53l l54l 156 1 had derived addition theorems for Slater-type 
functions by expanding them in terms of the complete and orthonormal functions ( 17.6b whose 
radial parts consist of generalized Laguerre polynomials. In the next step, Guseinov replaced 
the functions ( 17.6b according to ( 17.9b by Slater-type functions with integral principal quan- 
tum numbers, which are complete, but not orthogonal, and he also rearranged the order of 
the resulting summations. In this way, Guseinov essentially replaced expansions in terms of 
generalized Laguerre polynomials by power series expansions. Thus, the results of this article 
can in principle be used to check the legitimacy of Guseinov's manipulations. 

Since, however, one-range addition theorems for Slater-type functions are fairly compli- 
cated mathematical objects, the necessary determination of the decay rates and the sign pat- 
terns of the coefficients of the generalized Laguerre polynomials is very difficult. Fortunately, 
one can gain at least some insight by analyzing not the complicated addition theorems, but 
their much simpler one-center limits. In this way, it can for instance be shown that Guseinov's 
rearranged addition theorems for Slater-type functions with nonintegral principal quantum 
numbers do not exist for the whole argument range M 3 x M 3 . 
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A Divergent Series 

Divergent series have been a highly controversial topic that played a major role in the devel- 
opment of mathematical analysis (see for example the article by Burkhardt 1 26 1 or the very 
recent book by Ferraro ||44| ). and even now there is still a lot of active research on divergent 
series going on (see for example the recent monographs by Balser [4, 5|, Boutet de Monvel 
ifTTl . Candelpergher, Nosmas, and Pham [28|, and Sternin and Shatalov 1 123 1, or the review 
by Delabaere and Pham ll39ll ). 

As for instance discussed in articles by Barbeau |6|, Barbeau and Leah [7], Kozlov ll97ll . 
and Varadarajan 111 121 . already Euler had frequently used divergent series. Later, when the 
concept of convergence was better understood, Euler's admittedly somewhat informal treat- 
ment of divergent series was criticized, and a strong tendency emerged to ban divergent series 
completely from the realm of rigorous mathematics. This criticism culminated in Abel's fa- 
mous quotation from the year 1828, which very well expressed the prevalent attitude of most 
mathematicians during a large part of the nineteenth century and which can for instance be 
found in Littlewood's preface of Hardy's posthumously published classic on divergent series 
CD: 

Divergent series are the invention of the devil, and it is shameful to base on them 
any demonstration whatsoever. 

Although this disdain of divergent series was a temporary phenomenon of the nineteenth 
century, its consequences are nevertheless felt today. As emphasized by Suslov 11241 p. 1191], 
the standard university curricula in mathematical analysis were formulated in the middle of 
the nineteenth century. But this was a time when divergent series were wrongly considered 
to be essentially an aberration of the pre-rigorous mathematical past. As a consequence, 
divergent series and their summation are not part of the regular training of mathematicians 
and theoretical physicists, which in my opinion is totally inappropriate. In this context, a 
short article by Rubel 111 101 may also be of interest. 

At the end of the nineteenth century it was clear that the attempts of mathematical or- 
thodoxy to reject divergent series as unfounded had failed. Firstly, divergent series turned 
out to be too useful to be abandoned. For example, many special functions possess so-called 
asymptotic expansions which normally are factorially divergent inverse power series. In spite 
of their divergence, suitably truncated asymptotic series can provide excellent approximations 
at least for sufficiently large arguments. 

Secondly, the work of mathematicians like 

• Poincare, whose work in astronomy inspired him to formulate a mathematically rigor- 
ous theory of asymptotic series that typically diverge factorially, 

• Borel, who showed that factorially divergent series can be summed via Laplace-type 
integral representations, 

• Pade, who introduced his celebrated rational approximants that are often able to sum 
divergent power series, 

• Stieltjes, who showed that certain divergent series can be identified with continued 
fractions, 

ultimately led in the later part of the nineteenth century to a theory which made it possible 
to use divergent series in a mathematically rigorous way. Even more important from a prac- 
tical point of view is that their work showed that divergent series can actually be used for 
computational purposes if they are combined with suitable summation techniques. 
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In physics, divergent series are now indispensable. It is generally accepted that perturba- 
tion theory is the most important systematic approximation procedure in theoretical physics. 
But already in 1952, Dyson [42] had argued that perturbation expansions in quantum elec- 
trodynamics must diverge factorially, and since the seminal work of Bender and Wu on the 
perturbation expansions of anharmonic oscillators OUSl it has been clear that quantum me- 
chanical perturbation theory produces almost by default factorially divergent perturbation ex- 
pansions. A good source on divergent perturbation expansions in quantum mechanics and in 
higher field theories is the book edited by Le Guillou and Zinn-Justin |99| where many of the 
relevant articles are reprinted. 

Even if we agree that divergent series are now indispensable in physics, we nevertheless 
cannot expect that rigorously minded mathematicians are necessarily satisfied with the way 
divergent series are typically used in physics. Quite instructive is the following remark by 
Haldane from the year 1941, which is quoted in a book by Korner 11961 p. 426]: 

Cambridge is full of mathematicians who have been so corrupted by quantum 
mechanics that they use series which are clearly divergent, and not even proved 
to be summable. 

I think that the criticism, which is implicit in Haldane's remark, cannot be dismissed lightly. 
Let us for instance assume that we want to express a physical quantity by a divergent series, 
which for example may be a factorially divergent Rayleigh-Schrodinger perturbation expan- 
sion. Since divergent series are a priori mathematically undefined, we have to show that it 
is indeed possible to associate something finite - the value of the physical quantity - to the 
divergent series in a mathematically meaningful way, or to put it differently, we have to show 
that the divergent series is summable to the correct finite result by an appropriate summation 
method. 

Unfortunately, it is often extremely difficult to prove this rigorously. In such a situation, 
physicists tend to rely on their intuition and do not bother to try to formulate difficult proofs. 
Quite often, such a pragmatic approach is remarkably successful, but it should also be clear 
that mathematicians do not necessarily like that. Intuition can be misleading. Therefore, 
occasional unpleasant surprises and even catastrophic failures cannot be ruled out. 

B Sequence Transformations 

The fact that appropriate summation methods make it possible to use divergent series for 
computational purposes raises the question which of the numerous known summation tech- 
niques are best suited for the numerical evaluation of the divergent series occurring in this 
article. Typically, we are confronted with alternating series whose terms grow in magnitude 
like a fixed power of the index. Based on my own practical experience, I propose to use so- 
called sequence transformations, which are purely numerical techniques to transform a slowly 
convergent or divergent sequence {s„}™ =0 to another sequence {s'„}™ =0 with hopefully better 
convergence properties. 

For those interested in the history of sequence transformations, I recommend a monograph 
by Brezinski ||2TI . which discusses earlier work starting from the 17th century until 1945, as 
well as two articles by Brezinski [22, 23], which emphasize more recent developments. 

The basic assumption of all sequence transformations is that the elements of a slowly 
convergent or divergent sequence {s„ }" =0 , which could be the partial sums s„ = Y!l=o fl/t of an 
infinite series, can for all indices n be partitioned into a (generalized) limit s and a remainder 
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or truncation error r n according to 

s„ = s + r„, n e N . (B.l) 

If the sequence {s„}^ =0 converges to s, the remainders r„ in (IB.U can be made negligible by 
increasing n as much as necessary. But many sequences converge so slowly that this is not 
feasible. Increasing the index n also does not help in the case of the divergent series. 

Alternatively, one can try to improve convergence or accomplish a summation by com- 
puting approximations to the remainders r„ which are then eliminated from the sequence 
elements s n . At least conceptually, this is what a sequence transformation tries to do. 

With the exception of a few practically more or less irrelevant model problems, sequence 
transformations can only eliminate approximations to the remainders. Thus, the elements 
of the transformed sequence {*J,}^_q are also be of the type of ( IB. Il l, which means that a 
transformed sequence element s' n can also be partitioned into the same (generalized) limit s 
and a transformed remainder according to 

s'„ = s + r' n , neN . (B.2) 

The transformed remainders {r' n }™ =0 are normally different from zero for all finite values 
of n. However, convergence is accelerated if the transformed remainders {^,}™ =0 vanish 
more rapidly than the original remainders {r„}~ =0 , and a divergent sequence is summed if the 
transformed remainders vanish at all as n — > °°. 

Before the invention of electronic computers, mainly linear sequence transformations 
were used, which compute the elements of the transformed sequence {s|,}" = o as weighted 
averages of the elements of the input sequence {s n }„ = o according to 

n 

S 'n = E VnkSk- (B.3) 

k=0 

The theoretical properties of these matrix transformations are now very well understood and 
discussed in books by Hardy [78 ], Knopp 1941 . Petersen Q05 1, Peyerimhoff [ 106 1, Zeller and 
Beekmann 1 1 66 1 , Powell and Shah, and Boos [ 11 1 . Their main appeal lies in the fact that based 
on the work of Toeplitz 1 129 1 some necessary and sufficient conditions for the weights /x nJ t in 
( IB. 3b could be formulated which guarantee that the application of such a matrix transforma- 
tion to a convergent sequence {s n }™_ yields a transformed sequence {s'„}^ =0 converging to 
the same limit s = Sc°. 

There are also some so-called Tauberian theorems which show rigorously that certain or- 
ders of for example Holder's or Cesaro's summation method are needed to sum a divergent 
series whose terms grow in magnitude like a fixed power of the index ll78l Theorem 39 on p. 
95 and Theorem 49 on p. 103]. Since the terms of the divergent inner jj. series considered in 
Sections |4]and|6] all grow in magnitude like a fixed power of the index, it looks like a natural 
idea to use either Holder's or Cesaro's summation method for the summation of these diver- 
gent series. Unfortunately, the situation is not so simple. Since we also want to accomplish 

something useful if we only know a comparatively small number of numerically determined 

(ex) 

Laguerre series coefficients A„ , we should focus our attention on those purely numerical 
summation techniques that promise to be particularly efficient. 

From a purely theoretical point of view, regularity is extremely desirable and greatly facili- 
tates the formulation of nice mathematical proofs, but from a practical point of view, regularity 
is a serious disadvantage. This probably sounds paradoxical. However, Wimp remarks in the 
preface of his book 11631 p. X] that the size of the domain of regularity of a transformation 
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and its efficiency seem to be inversely related. Accordingly, regular matrix transformations 
are in general at most moderately powerful, and the popularity of most linear transformations 
as computational tools has declined considerably in recent years. 

Nonlinear sequence transformations have largely complementary properties: They are 
nonregular, which means that the convergence of the transformed sequence is not guaranteed, 
let alone to the correct limit. In addition, their theoretical properties are far from being com- 
pletely understood. Thus, from a purely theoretical point of view, nonlinear sequence trans- 
formations have many disadvantages, but they often accomplish spectacular transformation 
results which are clearly beyond the reach of regular matrix transformations. Consequently, 
nonlinear transformations now clearly dominate practical applications and - albeit to a lesser 
extend - also theoretical work. Detailed treatments of their theoretical properties and long 
lists of successful applications can be found in monographs by Brezinski [18. 1 9 20 ] , Brezin- 
ski and Redivo Zaglia |24), Cuyt (34] ED, Cuyt and Wuytack [37], Delahaye J40|, Liem, Lii, 
and Shih OHD, Marchuk and Shaidurov [Kj!, Sidi flU], Walz ED, and Wimp [163), or in 
reviews by Caliceti, Meyer-Hermann, Ribeca, Surzhykov, and Jentschura [27], Homeier [81 1, 
and myself 11361 ). 

In spite of their undeniable usefulness, nonlinear sequence transformations do not neces- 
sarily get the attention they deserve, in particular in the more theoretically oriented mathe- 
matical literature. For example, a (very) condensed review of the classical linear summability 
methods associated with the names of Cesaro, Abel, and Riesz can be found in Zayed's rel- 
atively recent book [165 Chapter 1.11.1], but the more powerful and computationally more 
useful nonlinear sequence transformations are not mentioned at all. Apparently, many math- 
ematicians still prefer to work on the theoretically very nice, but computationally at most 
moderately powerful regular matrix transformations. 

Nevertheless, there are encouraging signs that the situation is changing for the better, and 
there are now several books by mathematicians that describe how nonlinear sequence trans- 
formations can be employed effectively as computational tools. For example, the most recent 
(third) edition of the book Numerical Recipes 11081 now also discusses nonlinear sequence 
transformations (for a discussion of the topics treated there and for further details, see 015311 ). 

I can also recommend a recent book by Bornemann, Laurie, Wagon, and Waldvogel ifTBI 
on extreme digit hunting in the context of some challenging problems of numerical analysis. 
For this extreme digit hunting, the authors also use sequence transformations, whose basic 
theory is described compactly in their Appendix A. This Appendix is too short to provide a 
reasonably complete and balanced presentation of sequence transformation, but I think that a 
novice can benefit considerably from reading it. I also like the extremely pragmatic approach 
of the authors of this book, which is very uncommon among mathematicians. Probably, this 
is due to the fact that the authors are not primarily interested in the mathematical theory of 
sequence transformations: They only wanted to apply sequence transformations as computa- 
tional tools in order to obtain more precise results at tolerable computational costs. 

Then, there is a very recent book by Gil, Segura, and Temme [48 1 on the evaluation of 
special functions. It discusses in addition to various other computational techniques also 
Pade approximants, continued fractions, and nonlinear sequence transformations which all 
facilitate the evaluation of (power) series representations for special functions. 

My own research shows that nonlinear sequence transformations can be extremely useful 
in a large variety of different contexts. I applied them successfully in such diverse fields as the 
evaluation of molecular multicenter integrals of exponentially decaying functions 1 52 , 821 11221 
[1581[16D, the evaluation of special functions and related objects [84l l85l [1361 [1371 [1391 [1401 
1 14311 147ll 1481 1 1551 1 1621 . the summation of strongly divergent quantum mechanical perturba- 
tion expansions E E] [32 [33] [83] [85] [^ the 
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prediction of unknown perturbation series coefficients [8 83, 86l 11451 [146L and the extrap- 
olation of quantum chemical crystal orbital and cluster electronic structure calculations for 
oligomers to their infinite chain limits of stereoregular ^Mfli/-onedimensional organic poly- 
mers E2 EH usa. 

In view of all these examples, it is probably justified to claim that anybody involved in 
computational work should have at least some basic knowledge about the power and also 
about the shortcomings and limitations of nonlinear sequence transformations. 

Pade approximants [m/n]f(z) can be viewed to be a special class of nonlinear sequence 
transformation since they convert the partial sums f n (z) = T!l=o YkZ k of a (formal) power series 
for some function f(z) to a doubly indexed sequence of rational functions. As documented 
by the long list of successful applications in the monograph by Baker and Graves-Morris [3|, 
Pade approximants are now almost routinely used in theoretical physics and in applied mathe- 
matics to overcome problems with slowly convergent or divergent power series. It is, however, 
not nearly so well known among nonspecialists that alternative sequence transformations can 
at least for certain computational problems be much more effective than Pade approximants. 

The probably best known sequence transformation is Wynn's epsilon algorithm H164L 
which is defined by the following nonlinear recursive scheme: 

=0, 4 n) =s n , neN , (B.4a) 

4% = + {n J_ w > M e No . (B.4b) 

£ k £ k 

in) 

The elements e lk with even subscripts provide approximations to the (generalized) limit s 

of the sequence {s«}" = o to be transformed, whereas the elements e^+i w i m °dd subscripts 
are only auxiliary quantities which diverge if the whole process converges. A compact FOR- 
TRAN 77 program for the epsilon algorithm as well as the underlying computational algo- 
rithm is described in lfl36l Section 4.3]. In lfT08l p. 213], a translation of this FORTRAN 77 
program to C can be found. 

If the elements of the input sequence {s„}™ = q are the partial sums f„{z) — T!l=oYkZ k of 

the (formal) power series for some function f(z), then the epsilon algorithm with = f„ (z) 
produces Pade approximants to /(z): 

4i = [n + k/k] f (z), k,neN . (B.5) 

But Wynn's epsilon algorithm is not restricted to input data that are the partial sums of a 
(formal) power series. Therefore, it is more general and more widely applicable than Pade 
approximants. Moreover, the epsilon algorithm can be generalized to cover for example vector 
sequences. A recent review can be found in fBTI . 

Since the epsilon algorithm can be used for the computation of Pade approximants, it is 
discussed in books on Pade approximants such as the one by Baker and Graves-Morris 0, but 
there is also an extensive literature dealing directly with it. On p. 120 of Wimps book 11631 
it is mentioned that over 50 articles on the epsilon algorithm were published by Wynn alone, 
and at least 30 articles by Brezinski. As a fairly complete source of references on the epsilon 
algorithm, Wimp recommends Brezinski's first book lfl8l . However, this book was published 
in 1977, and since then many more articles on the theory or on applications of Wynn's epsilon 
algorithm have appeared. Thus, any attempt of providing a reasonably complete bibliography 
would be beyond the scope of this article. 

In a convergence acceleration or summation process, it is usually a good idea to try to 
use the available information as effectively as possible (possible exceptions to this rule are 
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discussed in B147ID . Let us assume that a finite subset {io^i, ■ ■ ■ ,$m} of sequence elements is 

available. Then, those elements e!$ of the epsilon table produced by ( 1B.41 > should be chosen 
as approximations to the limit of the imput sequence that have the highest possible subscript 
or transformation order. Thus, if m £ No is even, m = 2jj,,l use as approximation to the limit 
of the input sequence the transformation 01361 Eq. (4.3-4)] 

{so,su---,s 2f i} -» S^j , (B.6) 
and if m £ No is odd, m = 2fl + 1, 1 use the transformation 01361 Eq. (4.3-5)] 

{ji,S2,...,J2m+i} -> 4p- ( B - 7 ) 

With the help of the notation |xj for the integral part of x, which is the largest integer v 
satisfying V < x, these two relationships can be combined into a single equation yielding 
ESI Eq. (4.3-6)]: 

/ \ (m— 2[m/2J) ,yj q\ 

\ s m-2[m/2\T S m-2[m/2\ + l >••• i s mj — * £ 2[m/2\ ' 

Wynn's epsilon algorithm is an example of a sequence transformation that uses as input 
data only the elements of the sequence to be transformed. However, in some cases structural 
information on the dependence of the sequence elements s„ on the index n is available. For ex- 
ample, it is well known that the truncation error of a convergent series with strictly alternating 
and monotonously decreasing terms is bounded in magnitude by the first term not included 
in the partial sum and that it possesses the same sign as this term (see for instance ll94l p. 
259]). The first term neglected is also the best simple estimate for the truncation error of 
a strictly alternating nonterminating and thus diverging hypergeometric series 2^b(CKj j8; — z) 
with a,/3,z > 11291 Theorem 5.12-5]. Such an information on the index dependence of the 
truncation errors should be extremely helpful in a convergence acceleration or summation 
process, but a sequence transformations like Wynn's epsilon algorithm cannot benefit from it. 

A convenient way of incorporating such an information into the transformation process 
consists in the use of remainder estimates {o„}~ =0 . Because of the additional information 
contained in the remainder estimates, sequence transformations of that kind are potentially 
very powerful as well as very versatile. 

The best-known example of such a sequence transformation is Levin's transformation 
1 100 1, which is generally considered to be a very powerful as well as very versatile sequence 
transformation (see for example [24 8 fl l 12011121 II 13611 14911 and references therein): 

* ,(k\ (p+n + j) k - 1 s„ +j 
^ n) (P,s n ,(On) = LI , k,n£N . (B.9) 

K / h\ I R _l_ n J_ ii 1 



T( ivY* yfl +B+ j> 1 



Here, /3 > is a shift parameter. The most obvious choice is j3 = 1, which is exclusively used 
in this article. 

The numerator and denominator sums of JSJJ. (/3 , s n , (O n ) can also be computed recursively 
r ifl36l Eq. (7.2-8) - (7.2-10)] or in [149, Eq. (3.11)]): 

4" 3 = «», neN , (B.lOa) 

,(») _ r(n+l) +« + *)*-' („) 

L k+\~ L k (p+n + k+l) k k ' ^' neN °' (B ' 10b) 
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By choosing in ( IB. lOat either u„ = s n / (£>,, or u n = \/ CO„, we obtain the numerator and denom- 
inator sums of Levin's transformation dB.91 >. 

As discussed in more detail in 013611149 1. Levin's transformation is based on the implicit 
assumption that the ratio [s„ — s]/co„ can be expressed as a power series in l/(« + J3). A 
different class of sequence transformations can be derived by assuming that [s„ — s]/co„ can 
be expressed as a so-called factorial series, yielding [136, Eq. (8.2-7)] 



t(-iy 



k\ Q3 + « + s, 



n+j 



j~Q \jj + k)k- 1 COn+j 

Si H >{P,s n ,a>k) = - k — — , k,neN . (B.ll) 

y i-iy( ^ (P+n+j)k-i i 

7=0 



J J Q3+n + %_i On+j 

As in the case of Levin's transformation. j3 > is a shift parameter, and again, only j3 = 1 

^ (n) (j3,i„,(^)from^ 



is considered in this article. Formally, we obtain ^ n \p,s n ,CO„) from ^ n \p,s n ,(0„) if we 



replace in ( 1B.91 > the powers (/3 +n + j) k 1 by Pochhammer symbols (j3 +« + j)jfe_i. 

(n) 

The numerator and denominator sums of (j3 , s n , (0,, ) can also be computed recursively 
r if!36l Eq. (7.2-8) - (8.3-7)] or in G49] Eq. (3.12)]): 



5^ = u„, neNo, (B.12a) 
,(«+l) (j3+» + fc-l)(/3+n + fc) c(n) 
(l3+n + 2k-l)(p + n + 2k) 

The initial values u n = s„/co n produce the numerators of the transformation ( IB. lib , and the 
initial values u„ = l/(O n yield the denominators. 

Both Levin's sequence transformation ( IB. 91 ) as well as the related transformation (IB. 1 U 
utilize the information contained in explicit remainder estimates {©,,}~ =0 which should be 
chosen in such a way that the ratio [s„ — s]/co„ becomes a smooth function of n that can be 
annihilated effectively by weighted finite difference operators (see for example 11491 Sections 
II and IV] and references therein). Accordingly, the choice of the remainder estimates is 
of utmost importance for the success or failure of a convergence acceleration or summation 
process involving Levin-type transformations. In this respect, it may be interesting to note 
that a symbolic approach for the construction of asymptotic estimates to the truncation errors 
of series representations for special functions was recently developed in II 1 5 01 . 

On the basis of purely heuristic arguments Levin 111 0011 had suggested some simple re- 
mainder estimates which according to experience work remarkably well in a large variety of 
cases and which can also be used in the case of the sequence transformation ( IB. lib . But the 
best simple estimate for the truncation error of a strictly alternating convergent series is the 
first term not included in the partial sum l94l p. 259]. Moreover, the first term neglected is 
also an estimate of the truncation error of a divergent hypergeometric series 2Fo(a,b, —£) with 
a,b,z > ll29l Theorem 5.12-5]. Accordingly, Smith and Ford 1 120| proposed for sequences 
of partial sums of alternating series the following remainder estimate: 

(On = As„. (B.13) 

The use of this remainder estimate in [ 1 00 1 and dB.lU yields the following variants of the 
sequence transformations Q3 , s n , On) and ^ (j3 ,s n ,03 n ): 

4'Xp,s n ) = ^ n \p )Sn ,As n ), (B.14) 
5, (,,) (J3,.„) = y(" ] (l3,s„,As n ). (B.15) 
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If the elements of the input sequence {s„}™ =0 are the partial sums f„(z) = YJl=oYkZ k of the 

(formal) power series for some function f(z), thend^ (j3,/ n (z)) and 8% (j8,/„(z)) are ratios 
of two polynomials in z of degrees A: + n and fc, respectively (for a detailed discussion of these 
rational approximants, see 11491 Section VI]). 

The transformations d^\p,s n ) and in particular also Sj; n '(f},s n ) were found to be re- 
markably powerful summation techniques for divergent series with strictly alternating terms. 
Numerous successful applications of Levin's sequence transformation jB.9\ and of the re- 
lated transformation (IB- lib and their variants in convergence acceleration and summation 
processes are described in II 1491 pp. 1210 and 1225]. It may be interesting to note that re- 
cently 8^ n \p,s n ) has been used quite a lot in optics lfT2l 1X31 IT41 IT?! . 

In the case of the transformations (IB. 141 > and ( IB. 15b . the approximation to the limit with 
the highest transformation order is given by 

{s ,sx,...,s m+1 } -> Z$(P,sq), (B.16) 

where s["'(j3,s„) stands for either d^ l \p,s n ) or 8^ n '(P,s n ). 
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